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Abstract 

Digital  archiving  of  bioacoustic  data  provides  both  curatorial  and  scientific 
benefits.  To  realize  these  benefits,  key  system  requirements  must  be  satis¬ 
fied.  This  report  discusses  these  requirements,  and  describes  the  software 
tools  developed  by  the  WHOI  bioacoustic  laboratory  to  maintain  and  utilize 
an  archive  of  digitized  biological  sounds.  These  tools  are  written  in  standard 
C  code,  and  are  designed  to  run  on  PC-compatible  microcomputers.  Both 
the  usage  and  structure  of  these  programs  are  described  in  relation  to  the 
SOUND  database  of  marine  animal  sounds.  These  tools  include  software  for 
analog- to- digital  conversion,  text  header  maintenance,  data  verification  and 
interactive  spectrographic  review.  Source  code  listings  are  supplied. 
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1  Introduction 


Historically,  bioacoustics  laboratories  have  maintained  libraries  of  field  record¬ 
ings  and  prints  of  spectrographic  analyses.  Well-organized  catalogs  for  such 
collections  permit  routine  retrieval  of  tapes  or  spectrograms  in  the  inventory. 
With  computer  systems,  digital  technologies  can  be  incrementally  added  to 
these  collections.  A  more  comprehensive  review  of  potential  digital  appli¬ 
cations  suggests  that  fundamental  reorganizations  of  bioacoustic  programs 
are  appropriate. 

Bioacoustic  studies  are  increasingly  concerned  with  objective,  quantita¬ 
tive  analysis  of  sounds.  Such  results  can  be  difficult  to  obtain  using  spectro¬ 
graphic  prints.  Many  research  topics  require  analysis  of  thousands  of  sounds; 
this  would  be  virtually  impossible  without  dramatic  changes  in  techniques. 

Many  laboratories  are  now  confronting  the  difficulty  of  maintaining  or 
accurately  reproducing  aging  analog  tapes  in  their  archives.  With  analog 
recording  technologies,  gradual  degradation  of  the  data  cannot  be  avoided. 
Aging  magnetic  tape  media  and  tape  duplication  introduce  artifacts  and 
corrupt  the  original  recordings. 

Digital  technologies  provide  an  immediate  solution  to  the  archiving  prob¬ 
lem.  Entire  tapes  can  be  converted  to  digital  format  and  stored  on  a  variety 
of  media.  Although  no  medium  is  eternal,  some  optical  storage  technologies 
are  more  stable  than  magnetic  tape.  In  addition,  direct  digital  copies  can 
be  made  with  extremely  low  error  rates,  and  such  errors  can  be  detected 
and  removed. 

Important  benefits  can  be  realized  even  if  the  entire  inventory  cannot  be 
digitized.  Sound  sequences  can  be  digitized  and  stored  as  they  are  analyzed 
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in  the  course  of  research.  Re-analysis  can  proceed  with  the  digitized  copy. 
This  minimizes  the  handling  of  analog  tapes,  reducing  the  risk  of  damage  and 
degradation  from  repeated  usage.  Also,  it  should  be  much  easier  to  retrieve 
an  existing  series  of  digitized  sounds  than  to  locate  the  tapes  and  queuf’ 
the  sounds  for  re-analysis.  Digitized  sound  sequences  also  provide  “voucher 
specimens”  of  the  exact  acoustic  data  used  to  reach  scientific  conclusions. 

The  accumulation  of  digitized  sound  sequences  can  provide  new  opportu¬ 
nities,  by  increasing  the  scope  of  research.  For  example,  digitized  sequences 
produced  in  the  course  of  studying  individual  or  species-specific  repertoires 
can  be  combined  later  to  examine  the  structure  and  function  of  these  sounds 
in  broader  systematic  or  geographic  contexts.  To  realize  these  opportuni¬ 
ties,  a  laboratory  must  coordinate  computer  resources  and  solve  a  variety 
of  data  management  problems.  The  value  of  thousands  of  digitized  sounds 
is  dramatically  increased  by  a  flexible  means  of  organizing,  selecting,  and 
retrieving  them  based  on  associated  biological  and  environmental  data  (e.  g. 
species,  geographic  location,  season,  social  context,  individual  identity,  sex 
and  age).  We  have  achieved  this  capability  with  the  SOUND  databases  and 
associated  software  tools  (Watkins,  Fristrup  and  Daher  1991).  This  report 
describes  the  philosophy  and  implementation  of  our  developing  a  tool  set 
for  bioacoustic  research. 
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2  System  Overview 


Our  basic  resource  is  a  large  analog  tape  library  accumulated  over  more  than 
40  years.  The  reorganization  of  these  recordings  began  with  the  development 
of  a  database  to  improve  access  to  the  materials.  Each  tape  is  now  described 
in  an  entry  in  the  SOUND  database  (Watkins,  Fristrup  and  Daher  1991), 
which  includes  the  geographic  location,  time,  date,  species  present,  social 
context,  behavioral  observations,  and  other  relevant  information.  SOUND 
is  currently  maintained  with  INMAGIC  software  (Cambridge,  MA).  The 
structure  of  SOUND  makes  it  a  convenient  starting  point  for  new  research 
projects.  Database  queries  based  on  specific  research  needs  can  provide 
an  immediate  assessment  of  the  available  resources,  as  well  as  easing  their 
retrieval.  We  can  evaluate  the  research  potential  of  our  library  for  particular 
historical  periods,  geographic  regions,  phylogenetic  groups,  types  of  sounds, 
etc. 

Given  the  capacity  to  select  tapes  of  interest,  the  next  need  is  a  means 
of  converting  the  acoustic  data  into  digital  form  suitable  for  computer  data 
processing.  Two  distinct  applications  can  be  identified.  An  investigator  may 
need  to  browse  tapes  to  extract  particular  sounds  of  interest.  This  requires 
a  system  that  is  continuously  digitizing,  and  that  can  be  interrupted  by  the 
operator  to  save  the  most  recent  portion  of  the  signal.  Alternatively,  an 
investigator  may  wish  to  convert  entire  tape  recordings  into  digital  form. 
This  latter  application  requires  software  that  manages  the  rapid  transfer 
of  digital  data  from  an  analog-to-digital  converter  to  a  large  mass  storage 
device. 

Before  digitized  acoustic  data  are  accumulated  and  software  tools  are 
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developed,  standard  file  formats  need  to  be  defined.  The  file  format  should 
incorporate  a  text  field  containing  enough  information  to  ensure  proper  iden¬ 
tification  even  if  external  references  are  lost.  The  name  and  storage  location 
of  a  file  can  be  used  to  identify  sound  cuts,  but  the  embedded  information 
is  more  secure.  Digitized  sound  sequences  -  or  sound  cuts  -  are  of  little  use 
without  a  convenient  means  of  searching  for  and  retrieving  data  of  interest. 
A  second  database  is  an  obvious  means  of  organizing  sound  cut  information, 
and  it  can  be  closely  related  to  the  database  catalog  of  tapes.  Our  sound 
cut  database  is  called  SOUNDC  (Watkins,  Fristrup  and  Daher  1991). 

Other  requirements  for  maintaining  collections  of  digital  sound  cuts  are 
screening  the  acoustic  data  for  errors  in  digitizing,  cross-checking  the  text 
data  embedded  in  each  sound  cut  file  against  that  file’s  record  in  the  sound 
cut  database,  and  automatic  routines  that  manage  the  transfer  of  validated 
data  to  archival  mass  storage.  These  are  tedious  and  time-consuming  oper¬ 
ations;  efficient  and  reliable  methods  are  critical. 

A  program  for  reviewing  stored  sound  cuts  and  manipulating  the  data 
for  further  analysis  is  required.  This  should  provide  visual  and  aural  pre¬ 
sentations  of  the  sounds  as  well  as  the  ability  to  make  simple  measurements 
of  signal  features.  Additionally,  the  program  should  be  capable  of  exporting 
all  or  portions  of  a  sound  cut  in  convenient  file  formats  for  other  analyses. 

A  coordinated  system  that  meets  all  of  these  needs  is  described  below. 
Each  function,  along  with  the  structure  and  operation  of  the  relevant  soft¬ 
ware  tools,  is  discussed. 
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3  Digital  File  Format 

It  is  impoitant  that  a  bioacoustic  collection  standardize  the  format  of  their 
digitized  sound  files.  To  assure  compatibility,  we  have  adopted  a  header  and 
file  format  based  on  the  capabilities  of  a  standard  commercial  analyzer,  the 
Kay  Elemetrics  5500  Digital  Signal  Processor  (Pine  Brook,  N.  J.).  Their 
“5500”  (KAY)  format  is  used.  It  has  a  header  of  512  bytes;  the  data  follow 
in  the  form  of  16  bit  binary  integers  stored  in  2’s  complement,  low  byte, 
high  byte  format. 

The  header  is  a  mix  of  binary  and  ASCII  data.  The  binary  data,  written 
by  the  Kay  digitizing  equipment,  record  sampling  rate,  number  of  bits  of 
precision,  and  size  of  the  sound  cut.  The  text  field  is  used  for  storage  of 
arbitrary  notes  and  information.  In  the  header,  we  record  a  portion  of  the 
record  from  the  SOUNDC  database  that  describes  the  data  to  unambigu¬ 
ously  identify  the  source  of  the  material  and  the  digitizing  process. 

In  the  following  description  of  the  KAY  file  header,  the  bytes  in  the 
header  are  numbered  from  1  to  512.  Not  all  fields  are  listed. 

25-26:  ASCII  ‘1’  ‘2’  --  number  of  bits/sample. 

27-  :  ASCII,  null-terminated  string  --  number  of  samples  in  the  file 

65-71:  ASCII,  "5500SDM  --  Kay  5500  "trademark". 

121-122:  binary  integer  --  exponent  (base  10)  for  sample  rate. 
123-124:  binary  integer  --  mantissa  for  sample  rate. 

151-512:  ASCII  text:  ee  insert  SOUNDC  record  information 

File  names  code  the  identity  of  the  sound  cut.  The  first  two  characters 
are  the  last  two  digits  of  the  year  the  recording  was  made.  The  next  three 


characters  are  a  tape  number  within  the  year;  tapes  are  usually  numbered 
sequentially  as  they  are  cataloged.  The  last  three  characters  are  the  number 
of  the  digitized  sound  cut.  Because  more  than  999  cuts  nay  be  extracted 
from  a  single  tape,  we  use  a  36-base  numbering  system  for  the  cut  number. 
This  provides  over  40000  unique  cut  names  per  tape.  The  digits  are  0-9, 
followed  by  A  -  Z.  Thus,  the  cut  after  00Z  is  010.  The  remainder  of  t^e  file 
name  (the  MS-DOS  file  “type”)  is  .KAY.  Thus.  8721101A.KAY  is  the  46^ 
cut  (#01  A)  from  the  211th  tape  recorded  in  1987.  The  first  eight  characters 
of  the  name  also  provide  the  unique  identifier  for  the  corresponding  database 
record  in  SOUNDC. 
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4  Analog-to-Digital  Conversion 

We  satisfy  the  digitizing  requirements  of  our  research  with  two  modes  of 
analog- to-digital  conversion.  In  the  interactive  mode,  data  are  collected 
continuously  into  a  ring  buffer  of  random  access  memory  (RAM),  with  the 
oldest  data  being  overwritten  by  fresh  samples  as  they  are  acquired.  This 
cycling  is  interrupted  by  the  operator  when  a  signal  of  interest  is  com¬ 
plete,  and  the  contents  of  the  ring  buffer  (the  digitized  signal)  are  saved  on 
a  mass-storage  device  as  a  data  file.  In  practice,  we  use  real-time  sound 
spectrographic  processors  for  interactive  data  collection:  the  combination  of 
aural  and  visual  review  of  taped  material  promotes  more  effective  signal  re¬ 
view  and  identification.  In  addition,  the  interactive  mode  typically  supports 
higher  sampling  rates  because  the  data  are  stored  in  system  RAM,  which 
has  much  faster  storage  times.  About  a  minute  of  sound  can  be  stored  at 
the  highest  sampling  rates  in  our  systems. 

In  the  batch  mode,  an  extended  portion  of  a  recording  is  continuously 
digitized  into  a  large  file,  usually  on  a  hard  disk.  These  long  data  files  are 
subsequently  edited  or  processed  automatically  to  extract  signals  of  interest. 
Batch  digitizing  permits  larger  volumes  of  data  to  be  acquired  with  reduced 
operator  involvement.  We  routinely  digitize  tens  of  minutes  of  continuous 
sounds  at  a  sampling  rate  of  80  kHz.  Automatic  techniques  for  detecting 
and  characterizing  signals  are  planned  to  further  expand  our  data  processing 
capacity  (Fristrup  and  Watkins  1992,  K.  Christian  pers.  comm.,  J.  Buck 
pers.  comm.). 

The  Kay  Elemetrics  5500  Digital  Signal  Processor  (Pine  Brook,  N.  J.)  is 
a  commercial  interactive  sound  spectrograph  used  for  displaying,  selecting. 
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and  digitizing  sound  cuts.  The  Kay  can  sample  and  display  one  or  two 
channels.  Its  ring  buffer  has  8  megabytes  of  memory,  which  represents  51.2 
seconds  of  data  for  a  single  channel  at  the  highest  sampling  rate  available 
(81960  Hz).  Optional  hardware  and  software  permit  portions  of  the  data 
buffer  to  be  downloaded  into  an  ISA  bus  computer  running  MS-DOS.  The 
file  management  software  supplied  by  Kay  places  binary  information  (sample 
rate,  sample  precision  (bits),  data  length)  in  a  header  that  precedes  the  data. 
The  operator  can  choose  to  insert  notes  into  this  header  prior  to  downloading 
the  signal  data. 

VOICE  is  a  sound  spectrograph  developed  at  the  Woods  Hole  Oceano¬ 
graphic  Institution  (Martin,  Catipovic,  Fristrup  and  Tvack  1990).  VOICE 
consists  of  an  ISA  bus  microcomputer  running  MS-DOS,  augmented  by  ana¬ 
log  interface  and  signal  processing  boards.  Its  ring  buffer  is  limited  to  320 
kilobytes  of  memory,  and  the  maximum  sampling  rate  for  gap-free  data  is 
about  30  kHz.  At  the  maximum  rate,  the  buffer  holds  the  most  recent  five 
seconds  of  signal.  Any  portion  of  the  ring  buffer  can  be  saved  as  a  signal 
file.  At  present,  VOICE  does  not  support  placing  text  in  the  header,  but 
this  function  is  supplied  by  other  utilities  described  in  Section  5. 

CSTRM  is  a  batch  digitizing  program  that  runs  on  an  AT-compatible 
microcomputer  augmented  with  analog  interface  hardware.  At  present,  it 
is  configured  to  work  with  the  Canetics  PC-DMA12  analog  interface  board. 
Key  features  of  the  PC-DMA12  are  DMA  circuitry  for  both  A/D  and  D/A, 
100  kHz  maximum  sampling  rate,  anti-alias  filtering  and  programmable  in¬ 
put  amplifiers.  CSTRM  streams  data  to  hard  disk  in  the  KAY  format  at 
sampling  rates  up  to  83  kHz. 
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4.1  CSTRM:  A  Batch  Digitizing  System 

By:  Kurt  Fristrup 

System  requirements:  MS-DOS,  ISA  bus,  Canetics  PC-DMA12,  or  an 
equivalent  A/D  board  (with  modifications  to  the  code),  a  hard  disk  system 
(fast  and  large  capacity  preferable). 

Usage:  cstnn  <dest>  <*  64kbyte  blocks>  <sample  rate  in  Hz> 

Before  executing  CSTRM,  the  operator  should  use  a  utility  that  eliminates 
disk  fragmentation.  This  will  minimize  seek  time  required  by  the  hard  disk, 
and  permit  maximum  throughput.  The  operator  also  must  ensure  that  there 
is  enough  room  for  the  samples  requested.  At  present,  CSTRM  does  not 
verify  that  enough  disk  space  is  available.  We  have  found  it  convenient  to 
dedicate  a  physical  disk  or  logical  partition  for  streamed  data  and  move  the 
data  to  another  storage  location  after  each  session  to  clear  the  partition  for 
subsequent  use. 

The  command  line  arguments  are  the  disk/directory  destination  (a  stan¬ 
dard  DOS  path  description),  the  number  of  64  kilobyte  blocks  (32  kilosam- 
ples),  and  the  sample  rate  in  Hz.  We  have  digitized  up  to  190  megabytes  of 
continuous  sounds  with  CSTRM. 

Code  Description. 

CSTRM  utilizes  subroutines  adapted  from  those  supplied  by  Canetics 
with  their  PC-DMA12  analog  interface  board.  We  extensively  modified 
the  function  managing  continuous  acquisition  and  transfer  to  disk.  This 
program  is  compiled  using  the  COMPACT  model  of  memory  management 
(NEAR  code  branches,  FAR  data  pointers)  to  specify  memory  allocation 
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functions  that  use  far  pointers.  This  is  the  only  program  we  have  devel¬ 
oped  that  does  not  use  the  SMALL  model  (NEAR  code  branches  and  data 
pointers). 

The  first  statement  in  main()  calls  memreq()  to  obtain  a  full  64kb  page 
of  memory,  to  serve  as  the  buffer  for  the  DMA  data  transfers.  This  function 
allocates  128kb  of  memory,  and  returns  a  pointer  to  the  first  byte  (off- 
set=0x0000)  on  the  full  page  of  memory  embedded  in  this  block.  This  is  a 
bit  wasteful  of  memory,  but  it  does  not  limit  the  program’s  function. 

The  block  of  conditionals  involving  argc  and  argv[  ]  access  and  test  the 
command  line  arguments.  The  assignment  statements  referencing  sets—*  ini¬ 
tialize  a  data  structure  used  by  the  routines  controlling  the  analog  interface 
board  (A1B).  The  call  to  init.board()  initializes  control  parameters  in  the 
AIB,  and  kilLtimer()  disables  the  AIB's  timer  device.  This  ensures  that  the 
board  is  quiescent  prior  to  initiating  data  acquisition. 

The  subroutine  kfrdfile.contQ  performs  the  data  acquisition.  The  initial 
section  of  the  subroutine  initializes  several  pointers  that  monitor  DMA  sta¬ 
tus,  shift  the  data  to  5500  format,  and  transfer  the  data  to  disk.  The  skeletal 
5500  header  is  created  and  written  before  data  acquisition  is  initiated. 

The  general  strategy  involves  dedicating  a  64  kilobyte  page  of  memory  as 
a  DMA  memory  buffer,  and  setting  the  DMA  controller  to  cyclically  transfer 
data  into  this  page.  Thus,  the  DMA  controller  resets  itself  automatically 
when  a  full  page  of  data  has  been  acquired,  and  resumes  transferring  at  the 
beginning  of  the  page.  The  program  follows  the  DMA  controller's  progress, 
transferring  segments  of  data  to  disk  when  they  are  complete.  When  the 
page  is  divided  into  two  segments,  this  strategy  simplifies  to  a  ping-pong 
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buffer  system.  We  use  8  segments  at  present,  which  leaves  7/8ths  of  the 
buffer  available  to  buffer  data  acquisition  against  transient  delays  during 
disk  writes. 

The  logic  of  the  acquisition/ write  loop  is  relatively  simple.  The  high 
byte  of  the  DMA  pointer  register  is  polled  repeatedly,  and  all  of  the  new 
data  (up  to  the  DMA  pointer)  are  shifted  right  four  bits  to  translate  from 
Canetics  integer  format  (most  significant  12  bits  used)  to  5500  integer  format 
(least  significant  12  bits  used).  If  the  DMA  pointer  has  passed  the  end  of 
the  current  data  segment,  then  these  data  are  written  to  disk,  and  the 
pointers  marking  the  current  segment  are  updated  to  the  next  segment  in 
the  buffer.  The  last  segment  in  the  buffer  is  treated  differently,  because  the 
DMA  pointer  register  rolls  back  to  zero  when  it  increments  past  the  last 
byte  on  the  page.  This  mandates  a  slightly  different  test  to  determine  when 
the  last  segment  is  full.  This  code  also  updates  the  page  count  variable,  and 
terminates  sampling  when  the  page  count  equals  the  requested  sample  size. 

Program  termination  is  signalled  by  a  brief  sequence  of  audible  tones 
generated  by  the  computer’s  speaker.  The  program  also  prints  a  reminder 
that  the  DOS  clock  may  need  to  be  reset,  as  the  timer  tick  interrupt  is 
disabled  during  data  acquisition. 

A  slightly  faster,  multichannel  version  of  this  program  (CSCRM)  has 
also  been  developed  for  acoustic  localization.  This  version  does  not  produce 
a  KAY  formatted  file:  data  are  streamed  to  disk  in  the  native  format  of  the 
(Canetics)  PC-DMA12.  This  uses  the  most  significant  12  bits  of  a  16  bit 
word  to  store  each  sample,  while  the  KAY  format  uses  the  least  significant 
12  bits. 
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For  multichannel  recordings,  the  analog  interface  board  interleaves  data 
from  the  different  channels  into  a  single  file.  Accordingly,  an  additional 
program  (DUMSPLIT)  was  developed  to  demultiplex  the  multichannel  in¬ 
formation  in  the  composite  file  and  to  produce  separate  Kay  files  for  each 
channel.  In  addition  to  separating  the  data,  DUMSPLIT  shifts  the  sampled 
values  4  bits  to  the  right  to  convert  from  Canetics  to  Kay  sample  format. 
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5  Text  Header  Maintenance 


The  text  portion  of  the  file  header  identifies  the  source  of  the  data  and 
specifies  the  digitizing  settings.  This  information  should  be  complete  enough 
to  replicate  the  sound  cut  by  re-digitizing  the  original  recording.  Headers 
need  to  be  added  to  files  made  with  digitizing  instruments  that  do  not 
place  text  into  the  header  when  digitized  sounds  are  saved.  Furthermore, 
changes  in  database  structure  or  digitizing  procedures  may  require  altering 
the  format  of  the  text  already  in  place  for  a  large  number  of  files.  Therefore, 
a  program  that  permits  interactive  modification  of  individual  headers  is 
needed,  along  with  a  second  program  that  processes  large  groups  of  headers 
automatically. 

HEADEDIT  provides  the  interactive,  file-by-file  review  and  replacement 
of  the  header’s  text.  MASSEDIT  updates  the  headers  of  many  files  au¬ 
tomatically,  using  the  information  from  a  formatted  text  file  that  can  be 
imported  to,  or  exported  from,  the  SOUNDC  database  (Watkins.  Fristrup 
and  Daher  1991,  pp.  31-32). 

5.1  HEADEDIT:  Interactive  Header  Manipulation 

Bv:  Kurt  Fristrup 

System  requirements:  MS-DOS 

Uaaga:  haadedit  <filanama> 

HEADEDIT  reads  the  header  on  the  file,  and  displays  three  binary  fields 
(sample  rate,  number  of  samples,  number  of  bits  per  sample)  and  all  of  the 
text  information.  If  the  two  letter  codes  for  fields  used  in  the  text  database 
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are  recognized  in  the  header,  then  the  text  is  formatted  to  display  each  field 
separately.  Otherwise,  the  text  is  displayed  unformatted.  If  the  information 
is  correct,  the  file  may  be  left  unchanged. 

If  the  text  portion  of  the  header  is  to  be  changed,  type  [Ctrl-nj.  A  prompt 
appears  for  the  new  text,  which  may  be  typed  or  “pasted”  in  using  a  pro¬ 
gram  such  as  SIDEKICK  PLUS  (Borland  International,  Scotts  Valley,  CA.). 
When  the  text  is  complete,  an  [Esc]  terminates  the  entry.  HEADEDIT  ig¬ 
nores  more  than  the  first  362  characters  entered  into  the  header. 

The  new  header  can  then  be  accepted  [Ctrl-y],  retyped  [Ctrl-n],  or  you 
can  exit  the  program  leaving  the  header  unchanged  [Esc].  Our  use  of  Ctrl- 
y  and  Ctrl-n  for  yes  or  no  stems  from  problems  with  pasting  text  using 
SIDEKICK  PLUS:  when  we  attempted  to  paste  too  much  text  into  the 
header,  a  surplus  character  was  sometimes  interpreted  as  a  reply  to  the 
prompt.  If  you  type  [Ctrl-y],  the  new  header  is  inserted  in  the  data  file,  and 
the  program  exits.  If  you  type  [Ctrl-n],  the  program  cycles  back  and  you 
are  prompted  for  a  new  header.  [Esc]  aborts  the  program  without  changing 
the  data  file. 

Code  Description. 

HEADEDIT  opens  the  file  named  in  the  command  line,  and  reads  the 
first  512  bytes  into  a  buffer.  Integer  pointers  extract  the  mantissa  and 
exponent  of  the  sampling  rate  from  binary  fields  in  the  header.  The  number 
of  bits  per  sample  and  the  number  of  data  samples  are  stored  as  ASCII  text 
and  are  extracted  using  character  pointers.  The  text  portion  of  the  header 
(starting  at  byte  150)  is  searched  for  each  of  the  two  letter  codes  used  in  the 
SOUNDC  database  (e.g.  RN).  A  pointer  is  assigned  to  reference  each  code 
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found  in  the  header.  If  the  RN  (Record  Number)  field  code  is  not  found, 
the  text  is  displayed  verbatim.  Otherwise,  each  field  associated  with  a  two 
letter  code  is  displayed  on  a  new  line. 

The  input  section  takes  text  (from  the  keyboard  or  a  text  block  from 
SIDEKICK  PLUS)  and  places  it  into  the  header  buffer.  The  character 
combination  that  marks  the  end  of  a  line  (carriage  return,  line  feed)  is 
translated  to  an  underscore.  This  method  marks  the  end  of  each  line,  yet 
allows  most  software  packages  to  treat  the  text  information  as  one  continuous 
string. 

The  program  terminates  by  saving  the  new  header  information,  [Ctrl-y] 
or  retaining  the  old  header,  [ESC]. 

5.2  MASSEDIT:  Batch  Substitution  of  File  Headers 

By:  Kurt  Fristrup 

System  requirements:  MS-DOS 

Usage:  massedit  <taxtf ilanaaa>  <diak:path\> 

With  text  and  data  stored  in  distinct  locations,  and  many  researchers  con¬ 
tributing  to  the  database,  verifying  the  correspondence  of  text  and  data 
is  a  significant  concern.  Additionally,  the  format  and  contents  of  the  text 
database  may  evolve  with  time,  such  that  the  information  in  the  header 
becomes  out  of  date.  To  address  these  concerns,  we  use  the  same  text  file 
to  update  both  SOUNDC  and  the  sound  cut  headers;  alternatively,  we  up¬ 
date  the  sound  cut  headers  using  a  text  file  written  from  SOUNDC.  Thus, 
MASSEDIT  accepts  the  SOUNDC  formatted  text  records  and  inserts  por¬ 
tions  of  each  record  into  the  corresponding  sound  cut  file. 
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To  prepare  for  MASSEDIT,  a  text  file  must  be  created  with  the  partic¬ 
ular  format  expected  by  the  SOUNDC  database  program.  In  this  format, 
fields  within  a  record  begin  on  a  new  line,  starting  with  a  two  letter  code 
that  identifies  the  name  of  the  field.  Records  are  separated  by  a  dollar  sign 
symbol  ($)  placed  on  a  line  by  itself. 

The  first  command  line  argument  to  MASSEDIT  is  the  SOUNDC  text 
file  name.  The  second,  optional  argument  allows  the  operator  to  specify 
the  storage  location  of  the  KAY  files.  If  no  location  is  specified,  the  default 
directory  is  used.  As  the  program  executes,  screen  messages  indicate  flawed 
text  records  or  failure  to  locate  a  sound  cut  file  corresponding  to  a  text 
record. 

Code  Description. 

MASSEDIT  begins  by  retrieving  the  text  file  name  from  the  command 
line  argument  list,  and  opens  the  file  to  prepare  for  reading  records.  MASSEDIT 
reads  the  text  for  a  single  record  into  a  buffer  (stopping  on  encountering  an 
end  of  record  dollar  sign  ($).  The  buffer  is  then  searched  for  the  letters  RN 
(Record  Number),  and  the  eight  character  number  is  read.  If  the  RN  field 
is  not  successfully  retrieved,  an  error  is  reported,  and  the  first  fifty  charac¬ 
ters  of  the  record  are  displayed.  The  program  proceeds  to  the  next  record 
without  taking  further  action. 

The  eight  character  RN  code,  with  the  letters  “.KAY”  appended,  corre¬ 
sponds  to  the  name  of  the  digital  sound  cut  file  referred  to  by  the  record. 
MASSEDIT  prepends  the  di«k:path  information  to  the  file  name,  and 
attempts  to  open  the  file.  If  MASSEDIT  doesn’t  find,  or  cannot  open  the 
file,  an  error  is  reported.  This  flags  the  presence  of  a  text  record  for  which 
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no  sound  cut  was  found  (perhaps  from  accidental  erasure  of  the  sound  cut). 

If  the  file  is  found,  MASSEDIT  inserts  as  much  of  the  SOUNDC  text 
as  the  header  will  accommodate  (362  characters).  End  of  line  characters 
are  translated  into  underscore  characters,  to  avoid  a  conflict  in  the  KAY  file 
management  software.  If  the  record  contains  fewer  than  362  characters,  the 
remainder  of  the  text  buffer  is  filled  with  spaces. 

The  program  terminates  when  the  last  record  has  been  read  and  it 
reaches  the  end  of  the  text  file. 
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6  Data  Verification 


Duplication  of  text  and  data  flies  can  involve  large  numbers  of  files  and 
enormous  volumes  of  data.  Often,  these  files  may  exist  on  different  machines, 
perhaps  even  machines  that  use  different  operating  systems.  To  verify  that 
all  copies  are  identical,  we  have  devised  an  automatic  procedure  based  on 
two  public  domain  utilities:  CRC  and  DIF.  Versions  of  these  programs  are 
available  fur  many  different  machines.  This  procedure  can  be  applied  to 
any  subset  of  the  database,  and  on  different  machines  in  widely  separated 
locations.  We  present  here  the  sequence  of  instructions  that  is  appropriate 
on  MS-DOS  machines;  a  similar  sequence  exists  for  moot  other  machines. 

for  Xc  in  (*.kay)  do  crc  Xc  >>  dirnaas.crc 
sort  <dirnams.crc  >dirnams. srt 

now  copy  the  .srt  files  to  a  common  machine,  and  execute 

dif  dimaael.srt  dimaa«2.srt 

The  initial  CRC  command  generates  a  single  line  in  DIRNAME.CRC  for 
each  filename  that  satisfies  the  wildcard  specification.  In  addition  to  the 
name  and  size  of  the  file,  a  32  bit  CRC  value  (a  common  error  checking 
code)  is  generated  for  each  file.  The  sort  command  ensures  that  all  file 
entries  are  in  alphabetical  order,  and  the  DIF  command  will  identify  and 
print  differences  between  the  sorted  files.  This  will  find  file  omissions  as  well 
as  virtually  all  forms  of  file  corruption. 

MASSEDIT  will  identify  text  records  that  have  no  corresponding  sound 
cut  file,  and  update  the  headers  of  all  files  that  do  match.  It  does  not 
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detect  sound  cut  files  without  a  corresponding  text  record.  It  also  does  not 
check  the  validity  of  the  digital  data.  To  address  these  needs  we  developed 
KAYCHECK. 

KAYCHECK  scans  the  data  to  ensure  that  the  signal  was  not  clipped 
during  digitizing.  KAYCHECK  also  verifies  the  presence  of  a  SOUNE  J 
formatted  text  field,  and  checks  some  of  the  text  values  for  consistency  with 
the  binary  data  in  the  header.  KAYCHECK  provides  these  checks  on  every 
KAY  file  in  a  specified  directory,  and  it  generates  three  MS-DOS  batch  files 
to  assist  in  subsequent  diagnosis  of  flawed  files  transfer  of  validated  files, 
and  removal  of  archival  files  from  the  working  mass  storage.  On  exiting, 
KAYCHECK  displays  the  total  number  of  files  processed  and  the  number 
of  files  with  flaws.  The  batch  files  automate  tedious,  time-consuming  tasks, 
and  remove  the  possibility  of  faulty  copy  or  delete  operations  due  to  errors 
in  manually  issuing  these  commands. 

6.1  KAYCHECK:  Header  and  Data  Verification 

By:  Kurt  Fristrup 

System  requirements:  MS-DOS. 

Usage:  kaychsck  <diskdir«ctory>  [c] 

The  logical  structure  of  KAYCHECK  is  relatively  simple.  Text  output  files 
(KCOPY.BAT,  KDELETE.BAT,  KDIAGNOS.BAT)  are  opened  to  receive 
commands  for  copying  files  that  passed  and  diagnosing  files  that  have  faults. 

Each  KAY  file  (located  by  Turbo  C’s  findnextQ  function)  is  opened  and 
the  header  is  loaded  into  a  buffer.  The  record  number  (RN)  code  is  compared 
with  the  file  name.  The  sample  rate  (SR)  field  is  compared  with  the  binary 
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coded  sample  rate.  The  cut  size  (CS)  field  is  compared  with  the  quotient  of 
the  binary  coded  data  length  divided  by  the  binary  sample  rate.  If  any  of 
these  comparisons  detect  an  inconsistency,  an  error  message  is  displayed  on 
screen  and  a  line  is  added  to  the  KDIAGNOS  file  that  will  run  HEADEDIT 
on  this  sound  cut  file  to  diagnose  the  problem. 

The  data  in  the  file  are  screened  for  signal  clipping  that  might  have  oc¬ 
curred  during  digitizing.  At  present,  clipping  is  indicated  when  a  binary 
data  value  falls  below  1  or  above  4094  (the  maximum  A/D  sample  number). 
If  no  evidence  for  clipping  is  found,  the  minimum  and  maximum  sample  val¬ 
ues  are  checked  to  ensure  that  a  sufficient  fraction  of  the  sampling  dynamic 
range  is  used.  If  less  than  one-eighth  of  the  dynamic  range  is  used,  an  un¬ 
derflow  message  is  reported,  suggesting  that  the  signal  should  be  redigitized. 
In  either  case,  a  call  to  SIG  (refer  below)  is  inserted  into  the  KDIAGNOS 
file,  and  a  message  is  displayed  on  screen. 

If  all  tests  are  passed,  KAYCHECK  creates  am  entry  in  the  KCOPY 
batch  file  that  will  copy  the  data  file  to  a  new  drive  and  directory  and 
perform  a  binary  file  comparison  of  the  copy  with  the  original.  KAY  CHECK 
also  places  an  entry  in  the  KDELETE  file  to  assist  in  the  safe  recovery  of 
hard  disk  space  after  copying. 

Code  Description. 

KAYCHECK  begins  by  forming  the  wildcard  file  name  using  the  first 
command  line  argument  (drive directory  )  and  appending  “.KAY™ .  The 
three  batch  files  (KCOPY,  KDELETE,  KDIAGNOS)  are  also  opened  and 
assigned  to  FILE  variables.  The  Turbo  C  findfirstO  function  is  used 
to  locate  the  first  file  matching  the  wild  card.  These  steps  complete  the 


23 


preparation  for  the  principal  program  loop  which  follows. 

The  first  statements  in  the  loop  increment  the  file  count  variable,  open 
the  KAY  file  and  determine  its  length.  The  header  material  -  the  first  512 
bytes  -  are  read  into  a  character  array.  The  header  material  is  checked  for 
consistency  by  the  function  hsadbadO,  which  returns  a  non-zero  code  if  an 
error  is  found.  The  non-zero  code  causes  an  appropriate  screen  message  to 
be  displayed  and  an  entry  is  added  to  the  KDIAGNOS  batch  file  to  run 
HEADEDIT  on  the  KAY  file. 

If  no  header  errors  are  found,  the  data  are  scanned  (datscanO )  for  clip¬ 
ping  or  underflow  (insufficient  gain  during  digitizing).  These  errors  are  also 
indicated  by  a  non-zero  code  value  returned  by  datscanO;  an  appropriate 
screen  message  is  displayed  and  a  SIG  entry  is  added  to  KDLAGN'OS. 

A  counter  is  used  to  keep  track  of  the  number  of  files  with  problems, 
and  when  the  program  terminates,  the  total  number  of  files  scanned  and  the 
number  with  problems  are  reported. 
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7  Sound  Cut  Review 


To  review  stored  sound  cut  files,  a  flexible  tool  for  browsing,  editing  and  ex¬ 
porting  digital  acoustic  data  is  required.  In  addition  to  providing  a  graphic 
representation  of  the  signal,  this  program  provides  elementary  signal  mea¬ 
surement  functions  and  the  capability  to  edit  and  export  acoustic  data.  The 
most  challenging  aspect  of  developing  this  utility  has  been  the  wide  range 
of  signal  durations  and  amplitudes  that  it  must  handle.  The  graphic  dis¬ 
play  should  provide  a  reasonable  presentation  regardless  of  the  particular 
structure  of  the  signal. 

We  have  developed  a  utility  that  can  be  used  as  an  independent  program 
(SIG)  or  in  conjunction  with  the  text  database  (R-SIG,  a  TSR  “pop-up”). 
R_SIG  allows  immediate  display  and  manipulation  of  the  digitized  signals 
from  within  the  SOUNDC  database  program  by  swapping  it  out  of  memory 
and  loading  SIG  into  memory.  Both  versions  provide  reasonable  spectro- 
graphic  displays  with  the  default  settings,  and  can  work  with  files  of  any 
size.  Selected  portions  of  signals  can  be  expanded  for  more  detailed  display, 
or  exported  to  ASCII  or  MATLAB  (The  Math  Works)  formatted  files  for 
additional  analysis.  With  appropriate  hardware,  selected  portions  of  the 
signal  can  also  be  played  back  for  aural  review. 

7.1  SIG  and  RJSIG:  Sonagram  and  Waveform  Display  Util¬ 
ities 

By:  Kurt  Fristrup  and  Terrance  Howald 

System  Requirements:  MS-DOS,  ISA  bus,  numeric  coprocessor,  EGA  or 
VGA  video  graphics  system,  Canetics  PC-DMA12,  or  (with  modifications 
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to  the  code)  an  equivalent  D/A  board  with  DMA  capability,  amplifier  and 
speakers. 

Usage:  sig  filename. KAY  [dr]  [aa]  [e/v] 

Usage:  tsrint 

r.sig  <drive:path\  > 

R.SIG  is  subsequently  invoked  by  hitting  Alt-Esc. 

SIG  is  a  program  that  produces  spectrogram  and  waveform  displays  of  acous¬ 
tic  data  files.  Additionally,  portions  of  the  data  can  be  selected  for  more 
detailed  display  or  exported  into  ASCII,  MATLAB,  or  KAY  format  files. 
R.SIG  is  a  terminate-stay-resident  (TSR)  program  designed  to  work  with 
the  SOUNDC  database.  When  R.SIG  is  loaded,  any  file  whose  record  is  be¬ 
ing  reviewed  in  SOUNDC  can  be  displayed  with  SIG  at  the  touch  of  a  key. 
R.SIG  loads  a  small  (10  Kbytes)  keyboard  monitor  into  the  main  memory. 
When  a  hot-key  sequence  is  detected,  R-SIG  saves  the  SOUNDC  environ¬ 
ment  to  extended  memory  or  hard  disk  and  loads  SIG  into  memory.  The 
KAY  file  name  is  taken  from  the  SOUNDC  screen  by  R_SIG,  prepended 
with  drive :path  ,  and  passed  to  SIG.  When  SIG  is  exited,  R.SIG  restores 
the  SOUNDC  environment. 

SIG  can  process  long  sequences  (we  have  used  it  on  80  Mb  of  continuous 
data).  When  SIG  is  invoked,  the  screen  clears  and  the  program  begins 
scanning  the  data  to  obtain  scaling  parameters  for  the  waveform  and  spectral 
displays.  Just  prior  to  the  actual  display  of  the  signal  spectrogram,  time 
and  frequency  lines  are  drawn,  forming  a  grid  on  the  screen.  These  labelled 
lines  provide  a  convenient  means  of  interpreting  the  spectrograph  during 
the  drawing  of  the  display,  which  can  take  a  minute  or  more  for  long  files. 
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The  grid  is  erased  as  the  signal  spectrogram  is  drawn.  Cursors  are  used  to 
pick  out  time  and  frequency  information  once  the  spectrogram  display  is 
complete. 

The  completed  waveform  and  spectrogram  displays  include  basic  infor¬ 
mation  regarding  the  data  file.  The  top  line  contains  the  name  of  the  file, 
the  number  of  samples,  the  minimum  and  the  maximum  A/D  conversion 
numbers,  the  sample  frequency,  the  dynamic  range,  analysis  attenuator  val¬ 
ues  (the  latter  two  scale  the  conversion  of  spectral  power  to  color  values) 
and  the  FFT  size.  At  this  stage,  the  commands  tabled  below  can  be  used  to 
position  cursors  or  alter  the  display  parameters.  Four  cursors  are  available 
to  bracket  features  in  the  waveform  and  spectrogram  display.  These  provide 
a  means  of  scaling  the  image  and  measuring  time  or  frequency  intervals. 
They  can  also  be  used  to  focus  on  particular  sound  sequences  by  restricting 
the  range  of  data  displayed.  This  “zoom"  function  is  implemented  as  a  stack 
with  five  levels,  so  successive  commands  can  be  used  to  pick  out  local  details 
and  subsequently  return  to  the  more  general  display. 

Other  commands  include  toggling  the  spectrogram  display  on/off,  tog¬ 
gling  noise  compensation  on/off,  adjusting  dynamic  range  and  analysis  at¬ 
tenuator,  saving  a  section  of  the  data  to  a  new  file,  listening  to  the  data 
between  the  time  cursors,  changing  the  FFT  size,  and  exiting  the  program. 
The  format  of  an  exported  data  file  from  SIG  is  determined  by  the  exten¬ 
sion  specified  for  the  name  (.TEX,  .MAT,  or  other).  TEX  files  are  saved 
in  ASCII,  MAT  files  are  saved  in  MATLAB  (MathWorks)  format,  and  any 
other  extension  results  in  a  KAY  format. 
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COMMANDS  FOR  SIG 

left,  right  arrow 

move  a  time  cursor  left,  right. 

up,  down  arrow 

move  a  frequency  cursor  up,  down. 

Ctrl  left,  right  arrow 

move  a  time  cursor  5x  speed. 

Home  <keypad> 

toggle  front/back  cursors. 

PgDn  <keypad> 

zoom  in  to  new  display  bound. 

PgUp  <keypad> 

pop  back  to  previous  display  bounds 

Alt-P 

listen  to  data  between  time  cursors 

Alt-Fl 

multiply  playback  rate  by  two. 

Alt-F2 

divide  playback  rata  by  two. 

Ctrl  PrtSc 

export  data  between  the  cursors 

Esc 

exit  the  program! 

Alt-F 

toggle  spectrogram  display! 

Alt-N 

toggle  noise  compensation! 

Alt-D 

restore  default  parameters! 

Alt-F5 

multiply  fft  size  by  two! 

Alt-F6 

divide  fft  size  by  two! 

Insert  <keypad> 

increase  dynamic  range  by  3db! 

Delete  <keypad> 

decrease  dynamic  range  by  3db! 

Ctrl-PageUp 

increase  analysis  attenuator  by  ldb 

Ctrl-PageDn 

decrease  analysis  attenuator  by  ldb 

!  --  these  commands  are 

effective  during  screen  draw 
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Code  Description. 

The  code  for  SIG  is  split  into  several  files,  which  are  loosely  organized 
by  function.  The  core  subroutines  are  in  COMSIG.C,  including  mainQ. 
All  routines  related  to  FFT  processing  are  in  FFTFUNCS.C.  All  simple 
routines  related  to  screen  display  are  in  PLOTUTLS.C.  WHTIMIO.C  con¬ 
tains  a  few  input/output  routines  related  to  data  access  and  error  reporting. 
NDRAVVVLN.C  and  DRAWHLN.C  are  graphics  subroutines  that  utilize 
inline  assembly  code  to  maximize  line  drawing  speed.  NDRAVVVLN.C  is 
roughly  ten  times  faster  than  the  Turbo  C  library  function. 

MainO  copies  the  first  command  line  argument  (the  KAY  file  name)  into 
a  global  string  variable  filename,  sets  global  video  constants  and  spectro¬ 
gram  scaling  parameters,  and  passes  control  to  pstestQ.  This  structure 
reflects  our  adaptation  of  the  previous  generation  TSR  support  routines 
from  South  Mountain  Software  (South  Orange.  N.  J.).  This  required  that 
for  the  program  to  be  activated  as  a  TSR,  it  had  to  be  a  subroutine.  In  that 
version,  there  was  no  mainQ  in  COMSIG.C,  and  the  TSR  subroutine  used 
pstestO. 

PstestQ  switches  the  video  screen  to  graphics  mode,  opens  the  KAY 
file  to  read  the  data,  and  initializes  global  variables  in  the  FFTFUNCS.C 
file  using  functions  provided  for  that  purpose.  The  header  is  loaded  and 
the  sample  rate  is  decoded  from  two  integer  fields.  Next,  pstastO  calls 
gatac&leO,  which  selectively  scans  the  file  and  determines  appropriate 
scaling  factors  for  the  waveform  and  spectrogram  displays. 

GetscaleO  does  not  always  read  every  data  value  in  the  file.  Exceed¬ 
ingly  long  files  are  sampled  at  regular  intervals,  and  the  scaling  factors  are 
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estimated  from  these  samples.  GatacaleO  also  estimates  the  background 
noise  spectrum  for  the  sound  cut.  The  sampled  data  blocks  are  ranked  by 
acoustic  intensity,  and  those  that  fall  between  the  fifth  and  tenth  percentiles 
are  used  to  form  an  average  noise  power  spectrum.  This  spectrum  can  be 
used  to  adjust  the  spectrogram  display  such  that  the  background  noise  is 
“whitened”. 

After  scaling  factors  are  determined,  the  screen  display  commences. 
LineplotO  performs  the  waveform  and  spectrogram  display.  If  a  keystroke 
is  sensed  while  lineplotO  is  executing,  the  routine  aborts  and  returns  the 
keystroke  code.  This  permits  the  adjustment  of  the  scaling  factors  without 
having  to  wait  for  the  screen  display  to  complete. 

When  the  graphics  picture  is  complete,  a  brief  header  and  footer  are 
written  to  summarize  the  salient  characteristics  of  the  sound  cut.  Then, 
the  program  waits  for  a  keystroke  that  would  further  adjust  the  scaling 
parameters,  activate  time  or  frequency  cursors,  zoom  in  on  a  selected  portion 
of  signal,  or  export  a  portion  of  the  signal. 

Any  keystrokes  that  change  the  screen  display  (changing  dynamic  range 
or  analysis  attenuator,  modifying  FFT  size,  etc.)  cause  getcursesO  to  exit 
and  pass  control  back  to  the  principal  program  loop.  The  exit  code  is  used 
to  indicate  which  function  has  been  requested.  After  the  proper  parameters 
have  been  changed,  program  control  cycles  back  to  lineplotO. 

The  organization  of  the  graphics  display  is  constrained  by  the  specifi¬ 
cations  of  the  EGA  and  VGA  graphics  standards.  There  are  640  vertical 
lines  in  the  display,  so  the  portion  of  the  file  being  displayed  is  divided  into 
640  equal  segments.  If  there  are  fewer  than  640  segments  of  256  samples 
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each  in  the  file  (163840  samples),  the  envelope  display  is  completed  (no  sam¬ 
ples  skipped  during  min/max  search),  and  the  spectrogram  is  based  on  a 
256  sample  segment.  Because  the  spectrogram  buffer  includes  more  data 
than  each  waveform  segment,  the  spectrogram  will  tend  to  lead  the  wave¬ 
form  display.  In  an  extreme  display  of  one  sample  per  vertical  line,  features 
will  appear  in  the  spectrogram  display  nearly  one-half  screen  ahead  of  the 
waveform  display. 

If  there  are  more  than  163840  samples  in  the  file,  data  will  be  skipped 
in  computing  the  waveform  and  spectrogram  displays.  The  spectrogram 
display  becomes  the  average  of  several  subsegments  within  the  data  spanned 
by  one  vertical  line,  and  the  waveform  is  the  range  of  sample  values  observed 
in  those  subsegments. 

When  noise  compensation  is  enabled,  the  color  scaling  on  screen  is  re¬ 
lated  to  the  noise  floor  for  each  frequency  bin,  but  when  this  is  disabled, 
all  frequency  bins  are  scaled  using  the  same  standard.  The  noise  floor  is 
established  by  averaging  64  power  spectra  formed  from  segments  that  fall 
between  the  fifth  to  tenth  percentiles  in  observed  rms  values.  In  practice, 
noise  compensation  tends  to  bring  out  weak  signals  at  higher  frequencies. 

The  dynamic  range  parameter  adjusts  the  range  of  values  spanned  by 
the  fifteen  colors  used  in  the  display.  The  analysis  attenuator  parameter 
determines  the  range  of  values  placed  in  the  same  color  bin  as  the  loudest 
signal  found.  In  a  linear  relationship  between  intensity  (in  dB)  and  color 
value,  the  dynamic  range  sets  the  slope  of  that  relationship,  and  the  analysis 
attenuator  determines  the  offset. 

The  playback  feature  of  SIG  utilizes  the  D/A  capabilities  of  the  Canet- 
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ics  PC-DMA12  board,  or  an  equivalent  analog  interface  system.  The  initial 
playback  rate  is  set  equal  to  the  sample  rate.  The  user  may  increase  or 
decrease  the  playback  rate  from  a  minimum  of  640  Hz  up  to  a  maximum  of 
81920  Hz.  D/A  playback  is  initiated  by  entering  [Alt- P]  to  call  the  function 
start.dacO.  Parameters  passed  to  start.dacO  are  the  first  and  last  data 
points,  the  sampling  rate,  and  a  scaling  factor  for  the  sampling  rate.  This 
sound  playback  routine  executes  in  a  manner  similar  to  CSTRM,  with  the 
exception  that  start.dacO  performs  D/A  conversions  instead  of  A/D  con¬ 
versions.  Start.dacO  was  developed  from  routines  distributed  by  Canetics 
(Pasadena,  CA)  with  their  PC-DMA12  analog  interface  board. 

Start.dacO  allocates  a  64kb  page  of  memory  by  calling  memreqO  and 
sets  the  data  pointer  to  the  first  data  point.  The  program  then  determines 
the  number  of  samples,  the  user  selected  playback  rate  (sample  rate  times 
the  scaling  factor),  and  the  number  of  64kb  pages  to  be  converted.  Next, 
the  routine  initializes  the  AIB  and  executes  wtfile.contO  which  controls 
the  D/A  process. 

Wtfile.contO  examines  the  amount  of  sound  data  to  be  analyzed.  If 
this  amount  is  less  than  65535  (64kb),  then  the  data  are  read  from  the  hard 
disk  to  memory,  converted  from  the  KAY  integer  format  to  the  Canetics 
integer  format,  and  sent  via  DMA  to  the  AIB  in  16kb  blocks.  If  the  last 
16kb  block  is  not  complete,  the  remaining  portion  of  the  block  is  filled  with 
the  value  of  the  last  data  point.  The  PC-DMA12  operation  is  halted  after 
execution  of  the  last  block. 

If  the  amount  of  data  to  be  converted  is  greater  than  64kb,  then  the  first 
64kb  of  the  data  is  loaded  into  memory,  and  the  conversion  process  is  started. 
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After  the  memory  pointer  passes  the  first  16kb,  the  old  data  is  overwritten 
by  the  next  16kb  of  data  in  the  sound  sequence.  This  process  of  sending  a 
block  to  the  AIB  and  then  filling  that  block  with  new  data  continues  until 
less  than  16kb  of  remain.  The  last  block  is  handled  as  described  above,  and 
the  function  terminates 

The  code  for  R.SIG  consists  of  one  file.  Many  of  the  functions  used 
in  mainO  are  adapted  from  those  provided  by  South  Mountain  Software 
(Pasadena,  CA).  The  first  code  in  mainO  checks  to  see  if  the  switch  “rel" 
is  included  in  the  command  line.  If  the  switch  is  included,  the  program 
attempts  to  remove  R.SIG  from  memory.  The  program  reports  successful 
removal  or  failure  to  find  R.SIG  in  memory.  If  the  switch  is  not  present, 
the  program  checks  to  see  if  R.SIG  is  already  loaded  or  if  a  path  to  the 
data  files  is  included  in  the  command  line.  If  there  are  no  problems  with  the 
command  line,  swap.tsrO  decides  where  to  temporarily  store  the  SOUNDC 
data.  Then  inittsrO  loads  a  small  kernel  of  lOkb  into  the  main  memory. 

This  kernel  monitors  the  keyboard  for  a  hot-key  sequence.  When  the  key- 
sequence  [Alt-Esc]  is  detected,  runjigQ  is  executed.  This  routine  gets  the 
name  of  the  KAY  file  being  examined  in  SOUNDC  from  the  video  memory 
and  prepends  “SIG"  at  the  beginning  of  the  string  and  “.KAY"  at  the  end  of 
the  string.  SOUNDC  is  then  swapped  out  of  memory  and  SIG  filename. KAY 
is  executed.  Once  SIG  is  exited,  SOUNDC  is  swapped  back  into  memory 
and  the  user  can  continue  from  the  point  where  the  hot-key  sequence  was 
entered. 
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8  Conclusions 


The  system  that  has  been  developed  for  curating  and  analyzing  marine  an¬ 
imal  sounds  from  our  SOUND  databases  has  evolved  over  several  years. 
These  software  tools  have  satisfied  our  needs  remarkably  well,  and  we  an¬ 
ticipate  that  they  could  be  useful  to  others  for  similar  signal  database  man¬ 
agement.  Their  functions  are  summarized  as  follows: 

•  Maintenance  of  a  database  catalog  of  field  recording  resources  and 
their  contents. 

•  Integration  of  resources  for  interactive  and  batch  conversion  of  field 
recordings  to  digital  form  suitable  for  computer  data  processing. 

•  Adoption  of  a  standard  file  format  for  digitized  sound  sequences  that 
includes  a  text  field  for  identification. 

•  Maintenance  of  a  database  catalog  of  the  digital  sound  sequences  and 
their  contents. 

•  Creation  and  revision  of  text  information  attached  to  sound  data  files. 

•  Verification  of  the  presence  and  consistency  of  text  information  in 
sound  data  files. 

•  Validation  of  the  sound  sequence  data  (no  clipping,  sufficient  dynamic 
range). 

•  Coordination  of  means  for  finalizing  the  animal  sound  data,  copying 
validated  sound  files,  diagnosing  flawed  files,  and  removing  files  from 
working  mass  storage  once  they  have  been  archived. 
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•  Utilization  of  tools  for  reviewing  sound  data  files,  making  simple  mea¬ 
surements,  and  exporting  all  or  portions  of  files  in  convenient  formats 
for  additional  analysis. 

The  WHOI  marine  animal  SOUND  database  system  was  developed  be¬ 
fore  commercial  database  programs  were  available  to  store  and  manipulate 
arbitrary  binary  data  in  concert  with  text.  Several  programs  now'  offer  this 
capability,  including  the  capacity  to  develop  custom  routines  for  display¬ 
ing  and  manipulating  the  binary  data.  These  programs  would  remove  the 
necessity  of  coordinating  a  database  catalog  with  separate  archives  of  digi¬ 
tized  sound  sequences.  However,  they  would  not  offer  equivalent  latitude  in 
analytical  procedures.  A  critical  factor  in  evaluating  such  databases  is  the 
flexibility  allowed  for  custom  programming  and  the  ease  of  developing  such 
analytical  routines. 

The  evolution  of  research  requirements  emphasizes  the  need  for  flexi¬ 
bility.  The  SOUNDC  database  was  designed  to  encourage  selection  and 
retrieval  of  data  files  based  on  biological  criteria.  Our  recent  work  has 
focussed  on  automatic  characterization  of  marine  animal  sounds  (Fristrup 
and  Watkins  1992)  requiring  analysis  and  comparison  of  sounds  from  many 
species.  This  could  have  been  accomplished  by  performing  separate  select 
and  copy  operations  for  data  associated  with  each  species,  with  distinct  disk 
directories  dedicated  to  each  set  of  data.  However,  it  was  much  easier  to 
process  all  relevant  sounds  without  prior  sorting.  We  subsequently  linked 
all  of  the  numeric  results  to  their  respective  SOUNDC  entries,  and  used 
database  queries  to  sort  and  summarize  the  results. 
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We  discovered  that  INMAGIC  software  currently  used  by  the  SOUND 
databases  was  not  flexible  enough  to  conveniently  support  this  kind  of 
linkage  between  text  information  (which  remains  unchanged)  and  varying 
amounts  of  numeric  information.  Relational  database  programs  (such  as 
PARADOX,  Borland  International,  Scotts  Valley,  CA)  supply  this  flexibility 
by  providing  linkage  between  separate  tables  of  information.  The  SOUNDC 
data  was  translated,  therefore,  into  a  PARADOX  table  that  carried  the 
biological  information.  The  analytical  results  were  tabulated  separately, 
with  a  file  name  and  the  related  numerical  values  constituting  an  entry. 
The  SOUNDC  record  identifier  and  the  file  name  field  in  the  numerical  ta¬ 
bles  provided  the  means  of  linking  the  tables.  Categorical  summaries  were 
obtained  by  performing  queries  based  on  the  SOUNDC  information,  and 
reporting  linked  numerical  data  specified  by  PARADOX. 

It  is  not  possible  to  foresee  future  research  requirements,  but  extensive 
retooling  can  be  avoided  by  sufficient  scientific  insight  and  effective  data 
management.  The  time  invested  in  acquiring  and  annotating  the  digital 
bioacoustic  data  will  be  well  spent  if  the  resulting  system  facilitates  estab¬ 
lished  protocols  and  fosters  innovative  research. 
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9  Appendix:  Source  Code  Listings 


These  listings  represent  “working  code”;  they  are  not  necessarily  free  from 
bugs  or  structured  in  the  most  logical  fashion.  Continued  development  and 
modification  are  anticipated.  Most  of  the  code  for  these  programs  was  de¬ 
veloped  using  Turbo  C++  (Borland  International,  Scotts  Valley,  CA);  the 
exceptions  are  the  87FFT  library  (Microway,  Kingston,  MA)  and  some  video 
graphics  functions  that  were  coded  in  assembly  to  increase  screen  drawing 
rates.  We  used  standard  C  code  to  enhance  portability  (object-oriented  ex¬ 
tensions  were  not  used).  Compiler  switches  were  set  -  for  efficiency  -  to 
generate  80286  code,  use  a  numeric  coprocessor,  and  maximize  register  and 
jump  optimizations. 
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r 

Gutted  Canetics  routines  for  fast  streaming  to  disk 

MUST  BE  COMPILED  WITH  THE  COMPACT  MODEL!!!!!! 

kf  4/26/90 

*/ 

#include  <stdio.h> 

#include  <dos.h> 

#include  <stdlib.h> 

#include  < alloc. h> 

#include  < string. h> 

#include  <conio.h> 

#include  <fcntl.h> 

#include  <lib!80Ji> 

#  include  <stat.h> 

# define  CONTROLA  /*  add+lN  */  15876 

#  define  CONTROLB  /*  add+2*4  */  15880 

#define  OFFSET-ADD  /*  add+10*4  */ 15912 

#define  RESET-ADD  /•  add+ll*4  */  15916 
#define  DMAMODE-ADD  11 

#define  DMA-MASK-ADD  10 
#define  CLRPTR-ADD  12 

#define  AD-TIMER-MODE  122  /*  ad  is  timer  #1  in  this  mode  */ 

#define  DA-TIMER-MODE  52  /*  and  da  is  timer  #0  */ 

#define  TIMER-AD  15925 

#define  TIMER-DA  15924 

#define  TIMER-ADD  15924 

#define  TIMER-MODE-ADD  15927 

#define  STA  15904 

#define  STD  15908 

#  define  TOGGLE-DAC.ADDRESS  15888 

#define  BLOCK  SIZE  0x2000 

typedef  struct  settings  *psels; 
struct  settings  setsl; 
psets  sets; 

int  msb.  lsb,  array -size; 
char  c; 

unsigned  num-of-transfers, 
char  filename[30],  curfilname(30]; 
int  filel; 

int  fil«2 ; 

unsigned  runs  top.  inline,  gainbits,  inmode,  offset; 
char  fourloadflag,  flag; 
long  sets-tmp[llj; 
char  bits[9]  ; 

/* 

L'H-SEG  returns  the  segment  portion  of  an  anti-hugified  pointer 
UH.OFF  returns  the  segment  portion  of  an  anti-hugified  point®' 

*/ 

#define  UH-SEG(p)  ((  FP-SEG(  p  )  +  ((FP.OFF(  p  )  )  >>  4)  )  4  OxFOOO) 

#define  UH-OFF(p) 

((int)((((FP-SEG(p)+((FP-OFF(p))>  >4))fcOxOFFF)<<4)+(FP-OFF(p)&OxOOOF))) 
void  soundoff/ini  frq,  long  dur){ 
long  i; 
sound/  frq); 
for(i=0y<duru++); 

} 

int  far  *memreq(woid){ 
int  tempseg; 

void  far  ’death.to-nonbelievers  ; 
death-io-nonbelievers  =  farmalloc(0x20000L), 
if  (death-to-nonbe  lie  vers  ==  NULL){ 

pniX.fi" Corelrft  =  %ld  Need  %ld\n\r”,  coreleft(),  Ox20000L  ); 

putchar(l); 

exit/ 1 ); 
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} 

tempseg  =  UH-SEG(death-to-nonbelievers); 
if  (L'H-OFF(death-to.nonbelieva's)  !=  0) 
tempseg  +=  0x1000; 
return  (int  far  ")MK_FP(temp«eg,  0); 

} 

/* 

for  digital  I/O  and  muxes 

7 

#define  CONTROL-ADD  15876 
#define  OUTPORT-ADD  15932 
#defme  STROBE  15892 

/* 

gutted  liblSO.c  for  efficient  streamer 
kf  4/26/90 

V 

/* 

Thia  routine  initializes  the  PC'a  DMA  controller  prior  to  doing  conversions. 
It  masks  DMA  channel  1  on,  sets  the  channel’s  mode,  transfs-  direction, 
starting  address,  and  number  of  transfers  to  be  made. 

V 

void  init-dma_ad( unsigned  count,  unsigned  aeg){ 
outp(DMA_MASK-ADD.  0x05); 
out  p(  DMAMODE -ADD,  0x45); 
outp(CLRPTR-ADD,  0); 
outp(2.  0); 
outp(2,  0); 

outp(3,  count  &  OxOOFF); 
outp(3,  count  >>  8); 
outp(130,  seg); 
outp(  131,  seg); 
outp(DMA-MASK_ADD.  1); 

} 

/* 

This  routine  initializes  the  PC's  DMA  controller  prior  to  doing  conversions. 
It  masks  DMA  channel  3  on,  sets  the  channel's  mode,  transfer  direction, 
starting  address,  and  number  of  transfen  to  be  made. 

7 

void  init-dmajda( unsigned  count,  unsigned  seg){ 
outp(DMA-MASK-ADD,  0x07); 
outp(DMAMODE_ADD.  0x4B); 
outp(CLRPTR-ADD.  0); 
outp(6.  0); 
outp(6,  0); 

outp(7,  count  &  OxOOFF); 
outp<7.  count  >>  256); 
outp(130,  seg); 
outp<131.  seg); 
outp(DMA-MASK_ADD,  3); 

} 

void  initdma.ad.jont (unsigned  seg){ 
outp(DMA_MASK_ADD.  0x05); 
outp(DMAMODE-ADD,  0x55); 
outp(CLRPTR_ADD,  0); 
outp(2,  0); 
outp(2.  0); 
outp(3,  OxFF); 
outp<3,  OxFF); 
outp(130,  seg); 
outpf  131 ,  seg), 

outp(DMA-MASK_ADD,  0x01); 

} 

/* 

This  int  initializes  the  control  registers  on  the  PC-DMA  using  the 
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offset,  nil] nr  runstop,  datasise,  inmode,  gainbits,  and  filterstate 
global  variables  It  also  resets  the  board  before  and  alts'  to  insure 
that  the  board's  controlling  state  machine  is  in  the  proper  state  to 
start  conversions 

*/ 

void  init-board(psets  sets){ 
unsigned  temp; 

temp  =  sete->mline  +  (sets- > runstop  <<  5)  +  (sets->filterstate  <<  6)  ; 
outp( RESET-ADD,  0); 
outp(OFFSET-ADD,  sets->offset); 
outpfCONTROLA,  temp  +  (set*->datasize  <<  7)); 

outp(CONTROLB,  2  +  (sets->inmode  <<  2)  +  (sets->gainbits  <<  4)  +  128); 
outp(RESET-ADD,  0); 
if  (sets->datasize  =  =  0){ 

outp(CONTROLA,  temp  +  128); 
outp(CONTROLA,  temp  ); 

} 

if  (sets->dac2flag  ss  2) 

outp(TOGGLE-D  AC -ADDRESS.  0); 

}  /*  init-board()  */ 

/’ 

This  function  disables  all  the  counters  in  the  PC-DMA’s  on  board  timer 
device  Use  this  routine  before  changing  any  of  the  board's  operating 
parameters.  The  timer  activation  routine  should  be  the  last  thing  called 
before  conversions  start. 

*/ 

void  kill-timer(psets  sets){ 
outp(RESET-ADD.  0); 

outp(CONTROLB.  15  +  16*seta->gainbits  +  128); 
outpf  RESET-ADD.  0); 
if  (sets- >dac2flag  ==  2) 

outp(TOGGLE-DAC-ADDRESS,  0); 
outp(TlMER_MODE-ADD,  50). 
ouip(TIMER-MODE-ADD.  114); 
outpf  TIMER-MODE-ADD,  178); 

init-dma-da(0.  sets->seg); 
init.dma-ad(0,  sets->seg); 

} 

r 

This  routine  initializes  the  timer  device  on  the  PC-DMA,  setting  it 
up  to  generate  trigger  pulses  spaced  (spaces)  microseconds  apart. 

Note  that  the  A/D's  timer  is  operated  in  a  differait  mode  from  that 
of  the  D/A.  The  D/A  s  timer  is  a  straight  divide  by  n  counter  producing 
triggering  pulses  every  n  microseconds.  The  A/D's  timer,  on  the  other 
hand,  is  triggered  by  the  D/A's  timer  and  functions  as  a  programmable 
delay.  It  is  used  to  produce  a  triggoing  pulse  n  microseconds  after 
the  D/A  timer  produces  a  pulse.  This  was  int}ed  to  allow  Input  and 
Output  conversions  interleaved  at  a  rate  determined  by  the  D/A's  timer, 
and  spaced  (D/A  to  A/D  period)  by  the  A/D's  timer.  In  orda-  to  use  the 
A/D  alone  the  programmer  sets  up  the  timers  for  interleaving  but  sets  the 
D/A  to  be  triggered  by  software  control,  effectively  disabling  it  as  far 
as  the  timers  are  concerned. 

*/ 

r 

This  sets  up  the  timers  for  both  the  A/D  and  D/A  in  interleaved  mode, 
if  only  the  A/D  is  being  used,  set  the  inmode  variable  to  make  the  D/A 
controlled  by  software  strobe,  disabling  timer  control  of  it. 

V 

void  ini t-timer-ad( unsigned  spaces ){ 

outpf  TIMER-MO  DE-ADD.  DA-TIMER-MODE); 
outpf  TIMER-MODE-ADD,  AD-TIMER-MODE); 
outp<TIMER-DA,  spaces  4  OxOOFF); 
outp(TIMER-DA.  spaces  >>  8  )  ; 
outp(TIMER.AD,  2); 
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outp<TIMER-AD.  0); 

}  /*  init-timer-ad()  */ 

/* 

MU  up  the  program  for  streaming  to  disk 
modified  by  kf  4/26/90 

*/ 

void  kfrdfUe_cont(psett  sets,  citer  ‘curfilname )  { 

unsigned  jbuff,  curfil.  ipage,  blksz2,  blksiz,  blksizl,  pofset,  lastbuff; 

unsigned  far  “blkindx,  *irotptr,  “blkptr; 

kilLximer(sets); 

blksizl  =  (blk*z=BLOCKSIZE)  -  1;  blksz2  =  blksiz>>l; 

Iastbuff  =  (256(blksiz>>8))<  <8; 
curfil  =  openfcurfilname, 

O.WRONLY|O.CREAT|O.TRUNC|0-BLNARY, 
SJWRITE|SJREAD); 
blkptr  =  MK.FP(eet»->eeg<<12,0); 
memset  (blkptr, 0.51 2); 

•(blkptr  +  60)  =  2; 

•(blkptr  +  61)  =  10000/set»->waits; 

write(curfU, blkptr, 512);  /*  reserve  space  for  header  */ 

sets-  >itunode  =  l; 

init-board(sets); 

init-dma-ad-cont/tets-  >seg); 

pnntf(”reading  \n”); 

outp(0x21  JDxb9);  /*  kills  the  doe  clodt  interrupt  */ 

ipage  =  0;  jbuff  =  blksizl ;  blkindx  =  (irotptr  =  blkptr )-blksz 2; 

i  ni  t  -t  i  me  r  -ad  ( se  t  a-  >  w  aj  t  s ) ; 

outp(CLRPTR_ADD,  0); 

whtle(ipage  <  sets- > num-of-pag ) {  /*  ipage  is  64k  page  #  */ 

while( jbuff  <  0xffff){  /*  jbuff=ofs  lastbyte  in  current  buffer  */ 

while(inp(2)  jbuff  >=  (pofset  =  (un*igned)inp(2)<<8)) 
while(FP.OFF(irotptr)  <  pofset) 

*irotptr++  >>=  4; 
while(FP-OFF(irotptr)  <ss  jbuff) 

*irotptr++  >>=  4; 

/•  printT("%p  %p  %u\n"urotptr,  blkindx.  jbuff);*/ 

write(curfil,blkindx  +=blksz2.blksiz); 
jbuff  +=  blksiz; 

}  /*  first  n-1  buffers:  while  jbuff  <  Oxfflff  */ 
while(inp(2),lastbuff  <=  (pofset  =  (unaigned)inp(2)<<8)) 
while(FP-OFF(irotptr)  <  pofset) 

*irotptr++  >>=  4; 
while(FP.OFF(irotptr)  >  0) 

*irotptr+  +  >>=4; 

/*  printf(”9fp  %p  %u\n” .irotptr,  blkindx,  jbuff);*/ 

write(curfil.blkindx+  =  blk»z2  .blksiz) , 
ipage+  +  ;  jbuff  =  blksizl; 

} 

kilLtimer(sets); 

outp(0x21X)xb8);  /*  restores  the  doe  clock  interrupt  */ 
close/  curfil); 

}  /*  read-file xont{)  */ 
void  main(int  argc,  char  *argv(]){ 
int  far  *M-Ptr; 
long  sloop; 

M.Ptr  =  memreq(); 
sets  =  Jcsetsl, 
if  (argc  <  3){ 

print f(" usage:  cstrm  <filename>  <#  64  kb  block*>  <samplerate  Hz>\n"); 
exit(-l); 

} 

else  if  (argc  <  4) 

sets-  >  waits  =  20: 

else{ 
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if  (0.0==stof(srgv[3])){ 

printf(” problem  convming  %»  to  simple  rste\n",srgv[3]); 
exit(-l); 

} 

sets- >  wills  =  1 ,0e6/stof(irgv[3]); 
if  (14  >  sets- >  Wilts';  { 

printf("miximum  simpling  rite  for  DELL  316sx  system  is  T' 428Ht\n" ); 
exit(-l); 

} 

if  (stoi(srgv[2])) 

sets-  >num-of_peg  =  stoi(srgv(2]); 

else{ 

printf("problem  convoking  %s  to  number  of  64kb  blocks\n’',ergv[2]); 
exit(-l); 

} 

} 

printfC%d  microsecond  simpling\n%d  64kb  block»\n%e  minutes  of  simpling\n" , 
sets- >  waits  .sets- >num-of.pig, 

(douhle)65536.0e-6*aeta->num-of-pag’aeta->  waits/60.0); 
sels->dstisize  =  1;  /*  12  bit  iimpling  */ 

sets->mline  =  0;  /*  one  chsnnei  simpling  */ 

sets->runstop  =  0;  /*  no  multiplexer  cycling  */ 

sets->gsinbits  =  0;  /*  set  gsin  =  10  */ 

sets->inmode  =  3; 
sets-  >  offset  =  128; 
sets->filierstate  =  0; 
sets- >usingxmux  =  0; 
sets-  >dac2fiag  =  2: 

sets->seg  =  FP.SEG(  M-Ptr  )  >>  12  ; 

mun^of-transfers  =  sets-  >nuxrL-of_con  *  2; 

init-boird(sets); 

kilLtimer(sets); 

kfrdfile.cont(sets.  irgv(l]); 

for(sloop=400;sloop<1000;sloop= (sloop*  109)/ 100) 
soundoff(  sloop.25000) ; 

for  ( sloop=  1 000 ;sloop  >  400  ;sloop = ( sloop*  1 00)  / 1 09 ) 
soundoff(  sloop.25000) ; 
nosoundl ); 

printf("you  miy  need  to  reset  the  do#  clock!  use  ncc\n"); 

}  /*  miin  */ 
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r 

HEADEDIT:  inapect  and  edit  *.KAY  file  header* 

*/ 

#define  YES  0x19 
# define  NO  Oxe 
#include  <»tdio  h> 

#include  <atdiib  h> 

#include  <atnng  h> 

#  include  <conio  h> 
main(argc,argv) 
ini  argc; 
char  *argv[]; 

{ 

FILE  "fin; 
long  aamprate; 

char  ch,*hdbeg.  *hplr(40],  header(S13]; 
int  icmd,  i,  *imant,  *iexp; 

char  “atrcode  =  "RN  \0CU  \0NC  \0SR  \0CS  \0PL  \0SC  \0ID  \0AG  \0U  \0GS  \ 

\0GA  \0OD  \ont  \oda  \ogb  \ogc  \oot  \osh  \oos  \ona  \odi  \obh  \ohy  \orc  \org  \ 

\0RB  \0RL  \0ST  \0SL  \0SD  \0SP  \0FF  \0FS  \0SB  \0FE  \0AL'  \OLO 
header[51 2]=0; 
clr*CT< ); 

fin  —  fopen(argv(lj,"r+b" ); 
if  (argc  !=  2){ 

primR"  incorrect  command  tail\rai*e  head  edit  fname\n"); 

clr»cr( ); 

exit(2); 

} 

if  (fin  ==  N’L’LLK 

printf("%*  not  opened\n”,argv(l]); 
clr»cr( ); 
exit(  1 ). 

} 

fread(  header.  1 .5 12. fin); 

imant  =  (iexp  =  header  +  120)  +  1; 

aamprate  =  "imant; 

while(*iexp-) 

aamprate  =  aamprate*  10, 
icmd  =  0; 

hdbeg  =  headers  150; 

printf(”%a\n”  ,argv[l]); 

pnntf("aample  rate:  %ld\n",»amprate); 

printf("bita:  %c%c\n"  ,*(  header + 24).* (header +25)); 

pnnlf("*ample*  %a\n"  .beader+26); 

if(NLLL==*tratr(hdbeg,"RN  ")){ 

/*  dump  content*  on  acreen  for  diagnoaia  */ 
for(i=0;  i<512-150;  i  +  +  ) 
if(i%70) 

pnntfC%c",*(hdbeg+i)); 

elae 

printf("%c","(hdbeg+i)); 

}  /*  RN  not  found  */ 
el»e{  /*  RN  found  */ 
do{ 

/*  tag  all  recogruted  inmagic  code*  with  pointer**/ 
hptrliand]  ss  *tr*tr(hdbe;  atrrodel: 
if  ((hptr[icmdj  !=  hdbegltt(hptr[icmd)'=NULL)) 

*(hptr[icmd)-l)  =  0x7f; 
atrcode  +=  4; 
icmd+  +  ; 

)  w  hile  ( "atrcode  It  L  ( icmd  <  40 ) ) ; 
for(i=0;i<512-150u+  +  ) 

•witch  (*(hdbeg+i)){ 

/*  place  null*  before  all  inroagic  code*  found  */ 
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cu«  Ox7f  : 

*(hdbeg+i)=0, 

break; 

/*  convot  end  of  line*  to  underscore*  */ 
case  OxOa  : 
case  OxOd  : 

*(hdbeg+i)= 

break; 

} 

} 

/*  print  out  inmagic  field*  */ 
for(i=0;i<icmd;i+  +  ) 
if(NULL!=hptr[i]){ 

printf("  %*\n"  _hptr(i] ) ; 

} 

/*  if  change  desired,  read  from  stdin  and  write  into  the  file  */ 
printf("I»  the  header  correct  (ctrl-Y,  ctrl-N,  or  ESC)?"); 
do{ 

ch  =  toupper(Oxfll  getch( ) ) ; 
if(27==di) 

clr»cr().exit(0); 

}while(  YES!=ch  LL  NO!=ch); 
while)  NO==ch){ 

/*  ok,  we  want  to  change  the  header  information  */ 
clr»cr(); 

printf("pa*te  using  SideKick  now,  end  with  <E»c>\n"); 
for(i=0;i  <512-1 50. i+  +  ){ 

if(0xlb=  =  (ch=getch()i:0xff)) 
break; 

if((0xa==ch)[|(0xd==ch)) 

ch=V; 

putch(ch); 

*(hdbeg+i)=ch; 

} 

for(;i<512-150;*(hdbeg+i  +  +  )  =  '  '); 
clrscrf); 

for(i=0;i<512-150;i+  +  ){ 
putch(*(hdbeg+i)); 

} 

printf("U  thi*  correct’’  "); 
do{ 

ch=0xffi;getch(); 

if(27==ch) 

clr»cr(),exit(0); 

}while( YES!=ch  it  NO!=ch); 
if(NO==ch) 
continue; 

f*eek(  fin.  (long)  150,  SEEK  -SET); 
for(i=0u<512-150Jputc(*(hdbeg+i+  +  ),fin)); 

}  /*  while  NO  ==  ch  */ 
fcloaeffin): 

dr»or<); 

return(O); 

}  /*  main  */ 
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r 

MASSEDIT:  automatically  update  KAY  files  using  inmagic  TXT  file* 
for  every  RN  in  the  TXT  file,  look  for  a  corresponding  KAY  file 
and  insert  the  appropriate  header  information 

*/ 

#include  <stdio.h> 

#include  <stdiib.h> 

#include  <string  h> 

# define  TSIZE  512-150 
main(argc.argv) 
int  argc; 
char  "argvQ; 

{ 

FILE  "fin,  "fout; 
int  j,  i; 

char  ch,  headtext[TSIZE],  ofname[80]; 
if  (argc  <  2){ 

printf("  incorrect  command  tail\rmiae  maaaedit  fname\n"); 
exit(2); 

} 

if  (NULL  ==  (fin  =  fopen(argv[l],“r' ))){ 
printf("%s  not  opened\n"  ,argv[l]); 

exit(l); 

} 

do{ 

for(i=0,i<TSIZE  itit  !feof(fin);i  +  +  ) 
if(’\n'=  =  (ch=getc(fin))) 
headtext[i]  =  '.'; 
else  if('$'==ch) 
break; 

else 

headlext[i]=ch; 

if(TSIZE==i)  /*  clear  the  rest  of  the  record  */ 
while( !=getc(fin)&k'.feof(fin)), 

else 

for(j=ij<TSIZEieadtext[j-f+]  =  '  ); 

*hiler\n  !=getc(fin)  Lk  !feof(fm));  /*  clear  the  newline  character  */ 
strcpy(ofname.argv[2]); 

j=sscanf(headtext,’  RN  %[0-9A-Z*-i]“ ,ofaame+strlen(ofname)); 
if(»==j){ 

printf(" >%55.50s<  not  *canned\n"  headtext); 
continue; 

} 

strcat(ofname,"  kay"); 

/*  for(j=OJ<TSIZE;printf("%c"  ,headtext[j++])); 

printf(’<-\n");"/ 

if(NULL=  =  (fout=fopen(o£name,"r+" ))){ 
pnntf( "error  finding  %i\n',ohame); 
continue; 

> 

fseek(fout,(  long)  150.  SEEK-SET); 
fwrite(headtext,l  .TSIZE.fout); 
if(0!=fclo#e(fout)) 

prinlf('error  closing  %s\n"  ,ofnaroe;; 
memaet(headlext.O.TSIZE); 

}  while(!feof(fin)); 
fcloaef  fin); 
return(O). 

}  /*  main  */ 
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r 

KAYCHECK 

check  all  .KAY  file*  for  proper  annotation 
displays  each  file  that  fails  the  test 

7 

♦  define  HILIMIT  4095 

♦define  LOLIMIT  0  /*  limits  for  12  bit  integer  sampling  in  KAY  format  */ 

#define  HEADBYTES  512 

#define  NOTESTART  150 

♦define  RNNO  1 

♦define  SRNO  2 

♦define  CSNO  3 

♦define  RNBAD  -I 

♦define  SRBAD  -2 

♦define  CSBAD  -3 

♦define  SIZBAD  -4 

♦define  BUFFSIZE  4096 

♦define  BUFFSIZE2  8192 

♦  include  <stdio.h> 

♦  include  <io.h> 

♦  include  <dir  h> 

♦  include  <fcntl  h> 

♦  include  <time.h> 

♦  include  <math  h> 

♦  include  <string  h> 

unsigned  RN  =  ((uneigned)'N'<<8)|(unxigned)'R’; 
unsigned  SR  =  ((unsigned)'R'C  <8)|(unsigned)’S'; 
unsigned  CS  =  ((unsigned)'S'<<8)|(unaigned)'C; 

/* 

ERREXIT 

V 

void  errexitfchar  *msg){ 
fcloaeallf ); 

printfCerror:  %s\n",msg); 
getchf); 

exit(-l); 

} 

/• 

FOPENINP 

7 

long  fopeninpfchar  *fname.  int  *ifiie){ 

"ifile  =  open(fnaroe.O.RDONLY|O.BLNARY); 
while  (-1  ==  *ifile){ 

pnntf("Source  File  opening  did  not  succeed:  %s\n- .fname); 
printf(" Enter  file  name  for  input  file:  ”); 
if  (!scanf(”%s\n"/name)) 

errexit(" aborted  from  file  entry-); 

"ifile  =  open(£name,O.RDONLY|0-BINARY); 

> 

l*eek(*ifile,0,SEEK-£ND); 

return(teU("ifile));  /*  return  number  of  bytes  7 

} 

/’ 

DATSCAN 

reads  integer  data  from  binary  ifl.  checks  for  ove [range  errors 
computes  summary  statistics  and  writes  these  to  text  file  sS 

7 

ini  datscanfint  ifl.  FILE  *afi){ 

int  ii,  ibuflfBl'FFSIZE],  "iptr,  ired.  ival,  imin.  imax; 
imin=imax  =  (HILIMIT-LOLlMIT)/2; 
while(!eof(ifl)H/*  MAIN  DATA  LOOP  "/ 
ired  =  read(ifljbufTBL’FFSIZE2)>  >  1 ; 
for(ii=0.  iptr=ibufT;  ii  +  +  <  ired.  ){ 
if( LOLIMIT  >=  (ival="iptr+  +  )) 
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return(-l ); 

else  if  (HILIMIT  <=  ival) 
retum(-l); 
if(imin>ivzd) 
imin=ival; 
else  if  (imax<ival) 
imax=ival; 

} 

} 

if((imax-imin)"  (long  )16<(  long)  (HILIMIT-LOLLMIT)) 
retum(-2);  /*  underflow  */ 

else 

retum(O); 

} 

/* 

FINDLBL 

looks  for  a  two  character  sequence  in  a  null  terminated  string 

*/ 

char  "findlbl(char  "buff,  unaigned  kev){ 
unsigned  "iptr, 

while(iptr=buff.*'buff+  +  tt  *iptr!=key);  /*  move  int  ptr  bytewise  */ 
retum((*-buff?  buff  NULL)); 

} 

/* 

HEADBAD 

checks  for  inconsistent  recno.  sample  rate,  variant  cut  size 
returns  an  int  error  id 

V 

int  headbad(char  "fname,  char  "hdr.  long  fsiz){ 
double  cutsize; 

long  samprate.  hsrate,  nsamples; 
int  "imant.  "iexp; 
char  “keyptr,  sfieldfSO], 
imant  =  (iexp  =  hdr  +  120); 
samprate  =  ,  +  +imanl; 
sscanf(hdr+26."%ld"  .insamples); 
while(*iexp-) 

samprate  =  samprate*  10. 

if  (NULL’=(keyptr=fmdlbl(hdr+NOTESTART.RN))){ 
sacanf(keyptr."RN  %(0-9a-zA-Z]” .sfield); 
if  (stricmp(sfield. fname)) 

retum(RNBAD);  /*  not  equal  if  non  zero  */ 

} 

else 

return(RNNO); 

if  (NULL!=(keyptr=fijidlbl(hdr+NOTESTAKT.SR))){ 

sscanf( keyptr," SR  %lu*  ,t hsrate); 
if  ( hsrate  !=samprate) 

return/SRBAD);  /*  sample  rates  don't  agree  */ 

} 

else 

retum(SRNO); 

if  (NULL!=(keyptr=findlbl(  hdr  +  NOTEST  ART,CS))){ 
sscanf(keyptr,"CS  %le",t cutsize); 
cutwze-  =  (double)  nsamples /samprate; 
if(cuts;ze  >  0.1  ||  cutsize  <  -0  1)  /"  ern»  gtr  0  1  sec  "/ 
retum(CSBAD);  /"  cut  size  off  "/ 

} 

else 

return(CSNO); 

if  (nsamples  !=  (fsii-HEADBYTES)/2) 
return)  SIZB  AD); 
retum(0); 

}  /"  headbadf)  */ 
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main(int  argc,  char  *argv[]){ 

FILE  *fstat,*fcopy,*fdelete,*fcheck; 

char  header[HEADBYTES+l];  char  fnamerec(9),fname[80]; 
int  fin,  ifiles=0,  iwrong=0,  found,  headv&J; 
struct  ffblk  filerec; 
long  flength; 

strcpy(fname,argv[l]);  /•  must  terminate  with  a  backslash  */ 
atrcal(fname,"*.kay" ); 

if(NL'LL==(fcopy=fopen("kcopy  bat"  ,"w" ))) 
errexit(" error  opening  kcopy  file"); 
if  (NL'LL==(fdelete=fopen("  kdelete.bat"  ,"  w" ))) 
errexit(" error  opening  kdelete  file"); 
if(NULL=  =  (fcheck=fopen("kdiagnos  bat"," w"))) 
errexit(" error  opening  kdiagnos  file"); 
fprintf(fcheck,"ecfao  off\n"); 
found  =  !findfirst(£name,4:filerec,0); 
while(found){ 
ifilea+  +  ; 

header[HEADBYTES]=0;  /*  null  terminate  the  header  string  */ 

•trcpy (fnameArgvfl]);  /*  must  terminate  with  a  backslash  */ 

strcat(fname, filerec  ff-name); 

flength  =  fopeninp(fname,A:fin) 

lseek(fin.O.SEEK-SET); 

headval=read(fin,header.HE  AD  BYTES); 

fnamerec(8]=0, 

memcpy  ( fnamerec  .filerec .  ff  -name ,  8 ) ; 
if(headval=headbad(  fnamerec, header, flength )){ 

/*  pnntf("%s  %02d/%02d/%02d  %02d:%02d  ", 

filerec.ff.name, 

(filerec. ff-fd»tei;Oxlff)>  >5. 
filerec.ff-fdate&Oxlf, 

( filerec. ff-fdate  >  >9)+80. 
filerec.  ff  -ftime>>ll , 

( filerec  ff -ft imet0x3ff)  >>  5 ); ’/ 

iwrong+  +  ; 
switch  (headval){ 
caae  RNNO  . 

fpnnsf(fcheck."echo  RN  not  found\npause\n” ); 

break; 

caae  SRNO  : 

fprintf(fcheck  "echo  SR  not  found\npauae\n"): 
break. 

caae  CSNO  : 

fprintfffcheci.'echo  CS  not  found\npauee\n" ); 
break; 

caae  RNBAD  : 

fpn«f(fcheck."echo  RN  does  not  agree  with  filename\npause\n"); 
break; 

caae  SRBAD  : 

fpnr*f(fdieck,"echo  SR  does  not  agree  w'th  kay  value\npau»e\n" ); 
break; 

caae  CSBAD  : 

fprintf(fche<i,"echo  CS  off  by  more  than  0.1  sec\npauae\n" ); 
break; 

case  5IZBAD  : 

fpnntflfcheck  "ecno  file  length  number  of  samples \npau»e\n" )  . 

}  /*  switch  */ 

fprintf(fchedt,"headedit  %s\n"Aame); 
t  /“  if  headval  */ 

else  if  (headval=datacan(finjstat)){ 

/*  pnr*f("%s  %02d/(*02d/%02d  %02d  %02d  ", 

filerec.  ff  mame, 

(filerec  ff-fdatetOxlff)>> 5. 
filerec. ff-fdate&Ox  If. 
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(filer«c.ff-fdate>  >9) +80, 
fUerec.ff  Jtime?  >  1 1 , 

(filerec.ff  -ftime&:0x3ff)>>  5);*/ 
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if(-l  ==he*dval) 

fprmtf(fd5e<i,"echo  DATA  Ovemin\np«u»e\n” ); 

else 

fprintf(fche<i,"echo  DATA  Underflow \np»UM\n”); 
fprintf(fcheck,'’n*ig  %*\n" , fname); 
iwrong+  +  ; 

} 

el»e{ 

fprinlf(fcapy,"copy  %»  %%l\n”, fname); 
fprintf(fdelete,”del  %i\n” .fname); 

} 

cloee^fin); 

found  —  !findnext(fcfiler«c); 

}  /*  while  found  */ 

printf(”%d  file*  found,  %d  had  problexna\n’  dfiie*Jwrong); 

fcloaeall(); 

retum(O); 

}  /*  main  */ 
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#define  DRDEFAULT  27.0 
#define  DSDEFAULT  0.0 
#define  UNLOAD 
# define  SCRFILE 
#include  "whtimio.c” 

#include  "nsgbls.h" 

^include  ”d »ci" 

/"extern  unsigned  -stklen=0x2000;"/ 

/"  externals  used  to  permit  access  from  external  source  files  */ 
char  filename [N A MESIZEJ.  pathname[NAMESIZE], 

kolor«[VVSCALE],  header[HEADWORDS<<l); 

/*  globals  for  video  */ 

int  vidpage.  vscale=VVSCALE,  hscale=HSCALE,  bstgmode=  VBSTGMODE, 
botline=VBOTLINE,  ovidmode; 
int  idat[MAXFSIZE],  drsaturate; 
int  xfsize; 

Boat  scaler— 1; 
long  sampfreq; 
double  dynrange: 

/* 

global  float*  avgpow  assigned  only  in  getscale 

rdat  assigned.  referenced  many  place*,  colacale  assigned  in 

psteat  near  getscale,  maxpow  a  diagnostic  used  (potentially)  many  places 

■/ 

float  rdat[MAXFSIZE  +  2],  avgpow[(MAXFSIZE+2)>>l],  colscale,  maxpow; 

/* 

TICMARKS 

plots  ticmarks  to  aid  interpretation  during  drawing 

7 

void  ticmarks(lofig  he,  int  vmin,  int  fplot){ 
char  coir; 
double  tieval; 
int  i.tidine.tioexp; 

/* 

approx  5  horizontal  (time)  tics,  rounded  to  1  significant  digit 

V 

licexp=floor(  ticval=logl0((double)hc"HSCALE/(5“sampfreq))); 
ticval=floor(pow(  10  0,l>cval-ticexp  1+0.5); 
for(oolr=3,i=l ; 

HSCALE  >  (ticline= 

i"licval“»ampfreq*pow(10  0.(double)licexp)/(double)hc);i+ +  ){ 
draw  VLn(  ticli  tie,  vmjn.vscale.tcolr.O); 
gotcay((int  )((long)ticline"80/HSCALE).3); 
itoa(  (int )  tieval"  i.kolors.  10); 
wrilStringfkolors  .5,  vidpage); 

} 

gotaoty(40-8.2); 

wntStnngC seconds  xlO"  ,6, vidpage); 
gotoxy(43,l); 
itoa(ticexpjiolors,10); 
wmString(kolor»,6,vidpage); 
if  (fplot){ 

r 

approx  5  vertical  (freq)  tics;  assumes  never  <1  Hz  per  tic 

V 

ticexp=fioor(ticval=logl0(sampfreq/10  0)); 

ticval=floor(pow(10.0.ticval-ticexp)+0.5)*pow(10  0, ( double )ticexp); 
for(colr=3.i=l; 

256  >  (ticline=2*i*ticval"256/sampfreq- 1 )  J  +  +  ){  /*  xfsize  "/ 
drawHLn(  vscale-ticline.O.HSCALE- 1  ,colr); 
gotox  y  ( 70.  ( v  scale- 1  icline )  / 1 6 ) ; 
ltoat(long)ticval*  i.kolors.  10); 

•treat fltolors.”  Hz" ); 
writStnng(kolors.5,vidpage): 
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} 

} 

}  /*  ticmarks  */ 

/* 

LINEPLOTO 

divide*  the  time  domain  sequence  in  ifl  into  hinc  sized  cells 
starting  at  sample  ftrstd  (  ==byte  2*firstd);  hinc  precomputed 
and  draws  a  vertical  line  extending  from  the  maximum  value  to  the 
minimum  value,  the  vertical  scale  of  the  plot  remains  set  by  the 
maximum  found  in  the  file. 

*/ 

ini  linepiot(int  ifl,  long  firstd,  long  hinc. 

int  mnval,  ini  mxval.  int  fqplotjnt  xnoise) 

{ 

char  coir; 

int  vvmin,  vvmax,  vmax,  vmin,  vidmin,  vidmax,  i; 

long  akptween,  skpafter; 

int  ibuf,  nblkavg,  resid; 

int  hline.  noeof,  *iptr; 

float  vmult,  vadd.  “fptr, 

/*cls((char)0x70);*/ 

resid=xfrize;  /*  initialize  partial  screen  blank  marker  9-26-91  */ 
set  graphicsfbstgmode ) ; 

vidmin=16;  vidmax  =  (fqp]ot  ?  vscale  -  SPEKD1SP  :  vscale); 

ticmarks(hinc.vidmin.fqplot ); 

vmult= (float ((vidmin  -  vidmax)/(mxval-mnval); 

vadd  =  vidmax  -  vmult  ‘mnval; 

l*eek(ifl.((long)HEADWORDS+firstd)  <<  l.SEEK-SET); 
if(16  <  (nblkavg=hinc/256()  /*  xfsize  */ 
nblkavg=16; 
if(  nblkavg ){ 

skptween=(hinc-(\ong)256*nblkavg)/nblkavg;  /*  xfsize  “/ 
gain_set(gainjet(0.0)/(double)nblkavg):  /*  adjust  gain:  incoherent  */ 
} 

else{  /■  Inblkavg  •/ 
skptween=0; 

idatread(ifl,256.0ong)0.idat);  /*  xfsixe  "/ 

} 

skpafter=hjnc%(256-fskptween);  /*  xfsize  */ 

/* 

if  small  hinc  cause*  data  overlap,  get  initial  block 
subsequently,  the  old  data  is  shifted  right  and  new 
data  is  filled  in  on  the  left 

7 

for( noeof =  1 ,  hline=0,  hscale > hline ;  hline+  +  ){ 
vmax=-(vnun  =  MAXINT); 

fptr=rdat;  iptr=idat;  /*  used  in  both  halve*  of  if  hinc>xf*ize  */ 
if(nblkavg){ 

memaet(rdat,0,2  +  xfsize<  <2); 
for(ibuf=0;  ibuf+-f  <  nblkavg,) { 

noeof=idatread(ifLxfsizeskptween.idat); 

fminmax ( idat .  noeof,  fcvvmin.  tvvmax); 

vmin = min  (vmin.  vvmin);  vmax = max  (vmax.  vvmax); 

if  (fqplot) 

for(i=0;  i+  +  <  noeof;  *fptr  +  +  +  =  *iptr+  +  ); 
fptr=rdal.  iptr=idat;  /*  reset  pointers  */ 

}  /*  done  with  whole  blocks  */ 
if(  fqplot ){ 

flderoean<  rdat  .xfsize) ; 
makePspec  ( rdat  .xfsize ) ; 
iff  xnoise) 

normP  spec  ( rdat .  av  gpow ,  x  fsi  ze ) ; 
colrPipec(nial.kolor*,256.cal*cale).  /*  xfsize  */ 
drawVLnfhline. vscale- 256. vscale-1  Jtolors.l );  /*  xfsize  */ 
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} 

noeof=idatread(ifl,(int)skpafter.(!ong)0,idal); 
fminmax(idat,  noeof,  tvvmin,  icvvroax); 
vmin=min(  vmin,  vvmin);  vmax  =  max( vmu,  vvm»x); 

}  /*  if  hinc  >  xfsize  */ 
eu*{  /*  overlapping  segments  */ 
if  (resid>0){ 

fminmaxfidat,  hinc,  tvmin,  kvmu); 
if  (fqplot){ 

imean2fp(idat,rdat,  xfsize); 
makePspec(rdat,  xfsize); 
if(xnoise) 

normPspec(rdat,avgpow,256);  /*  xfsize  */ 
colrPspec(rdat,  kolors,256,colscale);  /*  xfsize  "/ 
drawVLn(hline,vscale-256,vscale-l  Jsolors.l);  /*  xfsize  "/ 
} 

memcpy(idat.idat+hinc,(256-hinc)<<l);  /*  xfsize  */ 
noeof=idatiead(ifl,hinc,(long)0,idat-f  256-hinc);  /*  xfsize  "/ 
resid-  =hinc-noeof; 

} 

else{ 

*kolors  =  0; 

draw  VLn(hlin«,l,v»calc,kolor8,0); 

} 

}  /*  if  hinc  >  xfsize..  else  */ 
if  (resid>0){ 

vmin=floor(0.5  +  vadd  +  vmult'vmin); 
vmax=:floor(0.5  +  vadd  +  vmult’vmax ); 
vn.  .i=min(vmin,vidmax); 
vmax = max  ( vmax ,  vi  dmi  n ) ; 

colr=0:drawVLn(hline.vidmin,vidmax,4:colr.0);  /*  colr=0  */ 
colr=12;drawV'Ln(hline.vmjn.vmax.i;colr.O).  /*  colr=12  */ 

} 

if(hioskey(l)){ 
if  (nblkavg) 

gain-sei(gain_set(0.0)*(double)nblkavg);  /*  reset  gain  */ 
retum(-l ); 

} 

}  /*  for  hline  */ 
relum(O); 

}  /"  function  lineplot  */ 

/* 

VCL'RSE 

*fln,  *Un  mark  the  begin  and  end  cursor  positions 

fdisp  and  ldisp  mark  the  begin  and  end  data  on  the  display 

*Un  and  ldjsp  mark  the  first  data  beyond  the  cursor/display 

imov  can  be  positive  or  negative,  and  it  has  been  multiplied  by  hacale  already 

*/ 

void  vcurse(long  *fin,  long  *lln,  long  fdisp.  long  ldisp. 

long  imov,  unsigned  lcurs) 

{ 

if  (lcurs){ 

if  (ldisp  <  (*lln  +=imov)) 

*Un=ldisp; 
else  if  (*Un  <  *ftn) 

*Un  =  *flru 

} 

else 

if  (fdisp  >  (“fin  +  =  imov)) 

*ftn=fdisp; 
else  if  (*Un  <  *fln) 

*fln=*Um 

} 

/* 
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HCURSE 

*fbar,  "lbar  mark  th«  begin  nod  end  frequency  ban 
imov  can  be  positive  or  negative 
upbar  ie  nonxero  for  upper  bar 

*/ 

void  hcune(int  *ub,  ini  *db,  ini  imov,  ini  upbar) 

{ 

if  (upbar ){ 

if  (SPEKDISP-1  <  (*ub  +=imov)) 

•ub=SPEKDISP-l; 
elaeif  (*ub<  *db) 

*ub=*db; 

} 

elae 

if  (2  >  (*db  +=imov)) 

*db=2; 

elae  if  ("ub  <  *db) 

*db=’ub; 

} 

/* 

TIMCAT 

timcat  adda  a  minutee/aeconda  aacii  airing  to  global  airing  kalon 

7 

void  timcai(long  jcurae){ 
ini  ijjt; 

char  una[Sj,  *cp, 

ltoa(jcurae/(60*aampfreq)  Jtolota+atrlen(kolon),10); 
atrcat(  kalon." 

ltoa((jcune/eampfreq)  %60,kolori-f  atriai  (kolora) ,  10) ; 
atrcat(kalon," ."); 

1  toa(  (j  curae*  ( long)  1 000 /aampfreq )%( fang )  1 000 ima .  1 0 ) ; 
if(0<(j=3-(k=atrlen(inie)))){ 

for(cp=ima+k-l.  i=0;  k>i+  +  ;  *(cp+j)=”ep,cp-); 
for(i=0,  cp=ima;  i+  +  <j;  *cp++='0’); 

"(im»+3)=0;  /*  terminate  with  null  */ 

} 

a  t  rcpy  ( kolora + a  t  rlen(  kolora )  ima ) ; 
at  real  (kalon,”  a"); 

} 

/* 

CRSWFUT 

aavea  repetitive  code  in  getcur*ea() 

V 

void  crawrit(long  i curae,  ini  roar,  int  la t curae,  ini  ip.  long  hnc){ 
ini  j,  a  min,  unax; 
gotcacy(20xow); 
at  rcpy  ( kolora,  “n=” ); 
ltoaficurae  jK>lora+atri«n(  kolora )  ,1 0) ; 
atrcaa(kolora,” ,  t=”); 
timcat  (icurae); 
atrcai(koiora," ,  a="); 
amajt=-(amin=MAXrNT); 

be«k(ip,(HEADWORDS+icurae)  <  <1 , SEEK-SET); 

j =nun(hncjdat  read  ( ip,  xfai»e,(  long  )0  jdat )); 
frninmax  ( idai  j  .it  amin  .Xcamaz ) ; 
i  ioafamin  .kolora + a  t  rien(  koto  re ) ,  1 0 ) ; 
atrcai(kolora,"  <  >“ ); 

itoaiamax  Jtolora+at  rim(kolora) ,  10) ; 
atrcat(kolora,"  pb="); 

ltoa(  aampfreq  *acal  er  holora+atrim(  kolora ),  10 ) ; 
atreat  (kolora*  Ht” ); 
if  (row) 

wntString(kolora,16+CURSCOL+latcurae,vidp^e); 

elae 
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tn  t  S  t  nng(  kolors  ,2  4  4-  C  U  RS  C  O  L-  Uterine ,  v  i  dpage ) ; 

}  /“  cr»wric()  “/ 

r 

GETCURSES() 

uiumei  that  the  screen  is  already  in  graphics  mode,  and  that  the  time 
domain  plot  has  already  been  drawn. 

defines  region  of  interest  using  begin  and  end  marltos  (vertical  lines); 
firstd  is  the  starting  data  sample,  Lastd  is  the  ending  data  sample 
hinc  is  the  horizontal  data  incronenl 

7 

int  getcuraes(long  “Mine,  long  “Ivine,  int  *ubar,  int  “dbar, 
long  ’fours,  long  *lcuis,  long  hinc, 
int  freqplot,  int  xnoise,  int  ifp) 

{ 

char  coir  ; 
int  iwait; 

unsigned  int  done,  ids,  j,  lastcuree,  last  bar, 
lastcurse=8;  /*  start  with  the  end  marker  */ 
lastbar=4; 

colr=CURSCOL+24; 

if  (hecalr  >  (j=(*lcurs  -  “fvline) /hinc)) 

drawVLn(j,16,vscale-l,fccolr,0);  /“  srrite  the  line  */ 
colr=CURSCOL+16; 
if  (0  <  (j=(*fcurs  -  *fvlme)/hinc)) 

drawVLn(j,16,vscale-l,tcolr,0);  /“  write  the  line  */ 
if  (freqplol){ 

colr=FCURSCOL+20; 
drawHLn(vscale-*ubar  ,0,539, coir); 
drawHLn(  v  scale- *ubar+ 1 ,0,639,colr); 
colr=FCURSCOL+16; 
drawHLn(vscale-*dbar  ,0,539, coir); 
draw  H  Ln(  v  scale-  “dbar + 1 .0 ,639, coir ) ; 

} 

if  (“fcurs  <  “fvline) 

“fcurs=“fvline; 
if  (“lours  >*lvline) 

*lcurs="lvline, 
done=done;  /*  aero  done  */ 

while(!done){  /*  loop;  active  curse  is  highlighted  as  LIGHTRED  */ 
wiule(0sdT  &  (ich=bioskey(0) ))  /“loop  until  kbd  sends  null  char*/ 

if(ich=  =  EXrTPROG)reium(EXnPROG); 
iwait=l;/“  borrowing  iwak  for  speed  shift  */ 
switch)  ids ){ 

case  0x7400  :  /*  ctrl-right  arrow  */ 

case  0x4d36  :  /*  shift-right  arrow,  move  marker  righs  5  “/ 
iwaits-iwait; 

case  0x7300  :  /*  ctrl-lrft  arrcsr  */ 

case  Ox4b34  :  /*  shift-left  arrow,  move  marker  left  5  */ 
iwait  *=8; 

case  Ox 4 bOO  ;  /*  left  arrow,  so  move  marker  left  “/ 
iwaits-iwait; 

case  0x4d00  :  /*  right  arrow,  so  move  marker  right  */ 
colr=CURSCOL+24; 

if  ((0  <  (j=((lastcurse’“lcurs:“fcvws)  -  *fvlkie)/hkic))l£i (hscak  >  j)) 
drawVLn(j,16,vscale-14iColr,0);  /*  xors  old  line  */ 
v  curse)  fours  Jcurs.  “fvline,  *lvline,(long)iwait*hinciastcurse); 
if  ((0  <  (j=((lasici*ee?*lcurs:“fci»s)  -  “fvline)/hkie))l£i£(hecale  >  j)) 
drawVLn(j,16.vscale-l<fecolr,0);  /*  xors  old  line  */ 
break; 

case  0x5000  :  /*  down  arrow,  move  bar  down  one  bin  */ 
iwaits-iwait; 

case  0x4800  :  /*  up  arrow,  move  bar  up  one  bin  */ 
if  ( freogslot )  { 

j=vscal«-(Ustbsr’’*ubar:*dbar). 
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drawHLn(j  +  +  , 0,639, FCURSCOL+20); 
drawHLn(j  ,0,639, FCURSCOL+20); 

hcurse(ubar,  dbar,  2*iw*it,  Leather), 
jsrvsc*le-(la*tbar?*ub*r:*dbar); 
drawHLjBU+  +  , 0,639, FCURSCOL+20); 
drawHLn(j  ,0,639, FCURSCOL+20); 

} 

break; 

case  0*4700  :  /*  keypad  del  key,  toggle  laatcroe  */ 
colr=16+(CURSCOL(CURSCOL+8)); 
if  ((0  <  (j=(*lcurs  -  *fvlme)/hinc))&fc(hsc*le  >  j ) ) 

dr*wVLn(j,16,vscale-l,fccolr,0);  /*  xors  old  line  */ 
if  ((0  <  (j=(*fcun  -  *fvline)/hinc))l£fc(h*cale  >  j)) 

drawVLn(j,16,v»cale-l  Acolr.O),  /*  ion  old  line  */ 
laatcune  ±8; 
if(frecg>lot){ 

colr=16+(FCL'RSCOL(FCURSCOL+4)); 
dr*wHLn(v*cal»*ub*r  ,0,639, coir); 
drswHLn(vscale-*ubar+l,0,639,colr); 
dr*wHLn(vscale-*dbar  ,0,639, coir); 
drawHLn(vscale-*dbar+ 1 ,0,639,  oolr); 
last  bar  ±4; 

} 

break; 

case  0x6800  ;  /*  Alt-Fl,  increase  freq  by  a  factor  of  2  */ 
scaler  *=  2; 

if  (sampfreq*ecaler  >  81920) 

{ 

printf(-\l7-); 
scaler  /=  2; 

} 

break; 

case  0x6900  :  /*  All-F2,  decrease  freq  by  a  factor  of  2  “/ 
scaler  /=  2; 

if  (sampfreq*scaler  <  640) 

{ 

pnnif(*\*r ); 
scaler  ’=  2; 

} 

break; 

case  0x1900  :  /*  ait-p,  playback  data  between  cursox  */ 
start  ulac  ( ifpjcun  irurs  aampfreq  scaler ) ; 
break; 

case  AAPLUS  : 
case  AAMINUS  : 
case  DRPLUS  :  /•  */ 
case  DRMINUS  :  /•  */ 

case  FDEFAULT  :  /•  alt-d,  restore  sonagram  default*  */ 
case  XNOISE  :  /*  ait- a.  toggle  noise  suppression  */ 
case  WRTFILE  ;  /*  save  data  within  cursors  in  new  file  */ 
case  PSHGSTACK  :  /*  control  page  down,  redisplay  with  new  limits  */ 
case  POPGSTACK  :  /•  control  page  up,  pop  previous  cursors  */ 
case  FREQDISP  :  /*  alt-f,  draw  sonagram  */ 
case  EXITPROG  :  /*  escape,  done  moving  markers  exit  */ 
case  FFTNCRSE  :  /•  alt-FS,  increase  fft-siae  by  a  factor  of  2  "/ 
case  FFTDCRSE  :  /*  alt-F6,  decrease  fft-sise  by  a  factor  of  2  */ 
done  —  1 ,  /•  terminates  outer  while  loop,  action  taken  */ 

/*  in  the  petest()  while  loop  */ 

}  /*  switch  ich  */ 
iwait=0: 
if  ( ’done) 

while)  'btcakey  ( 1 ) ) 

tf(5000<  =  iwail  +  +  ){ 
clr  line  (bot  line); 
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clrline(O); 

•trcpy  flsolor*  .filename ) ; 
writStringf  kolon.  5,  vidpage) ; 

cmwrit(*fcur*,0.1a*tcur»e ,ifp.hine);  /*  from /upper  cur»or  */ 
if  (xnoiae) 

writStringf"  noiae  refermced",  5,  vidpage); 
gotaayfO, hotline); 
if  ffreqpiot){ 

lto*((*ubar>  >  1  )**amp£req/SPEKDISP.kaiori,10); 
writString(kolom.FCURSCOL+16+laetbnr,vidpage); 
wri  t  CW  ( ’ ,  5 ,  v  i  dpag  e ) ;  euraor*  c  ( 1 ,  vi  dpage ) ; 
ltoa((*dbar>>l  J'tnmpfre^SPEKDISPJtalon.lO); 

wri  t  S  i  ring  ( kolors  ,F  CU  RS  C  0  L+  20-  Us  t  bar ,  vidpuge ) ; 
writStringf"  ="  , 5,  vidpage) ; 

ltoa(((*ubar>>  1  )-(*dhar>  >  1  ))%*mp£re<i/SPEKDlSP,kolor»,10); 
writString(kolor»,FCURSCOL+ 18,  vidpage); 
writStringf"  H*"  ,5,  vidpage); 

cm  wri  t  ( "lcur* ,  bo  l  line  4a»  t  eune  .ifp  June ) ; 

•trcpyfkolom," ,  dt="); 
timcat(*lcur»-*fcur»); 

writString(kolor»,20+CURSCOL,vidpage); 

break; 

}  /•  if  j++=iooo  v 

}  /*  while  (!done)  */ 

*f.line=*fcur*; 

*lvhne=*lcum, 
re  turn  fieri); 

}  /*  function  getcume  “/ 

r 

rWKTFILE 
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int  wrtfUefint  ifile,  long  fdiip,  long  Idiip)  { 
int  i,  nrd,  itext.  "iptr; 
long  iaamp.  *mh; 
char  fname[40],  *cptr; 
char  mathead[50j; 

FILE  *ofl; 

/•claf(char)0*7D);*/ 

•et  graphic*  (batgmode ) ; 
gotcwy(O.O); 

writStrmg("oulput  file  name:  " ,Oxf,  vidpage); 
atraetffname.O); 

while  (0=  =  *fname)  /*  read  String  return*  NULL  if  eacapcd  */ 
readS  t  ringf  fname  ,0 ) ; 

if  (itext=(NULL  !=atmtr(£name,"  lex" ))) 

ofl=fop>en(£name,'w" ); 

elae 

oil =f open  (fname.  "wb" ); 
if  (ofi==NULL){ 

writStringf "deatinatioo  File  opening  did  not  rucceed..  "  .OriX, vidpage); 
rettrn(-l); 

laeak(ifile,(H£ADWORDS  +  fdiap)<<l. SEEK -SET);  /•  aame  for  all  format*  */ 
if  (itcit){  /*  writing  in  aacii  */ 

forfiaampsfdiap;  iaamp<ldi*p;  iaamp-f  =nrd){ 

nrd=)datread(if3e.nrt±=min(xfai*e  Jchrp- iaamp).  (long  )Ojptr=idat); 

for(i=0;  i-f  +  <nrd;  ){ 

itoa(*iptr+  +  .ep4r=  fname, 10);cptr-; 
whJef  *+-f cptr)  /*  my  fwriteatring  function  */ 
putc((im)*cptr,ofl); 
putc('\n’.ofl); 

}  /*  for  i  */ 

}  /*  for  iaamp  */ 
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}  /*  if  itext  */ 

elae  if  (NULL  !=strstr(fname,"  mat" )){  /*  wriun*  matlnb  •/ 
mh=mathead, 

*mh+  +  =  40; 

*mh+  +  =  l; 

*mh+  +  =ldup-fdifp; 

*mh+  +  =0; 

ultaa(aampfreq  .k  alor* ,  1 0) ; 

*mh=alrlen(kolor»)  +  2; 

"(mathead+20)  =  V; 

•trcpy(mathead+21.kok>n); 
fwnte(mstbead.  1 ,21  +strien<  mathead+  20)  ,ofl); 
for(isamp=fdisp;  iaamp<ldisp;  isamp+=nrd){ 

nrd=idatread(if3e,nrd=min(xfsi»e>ldUp-isamp),(long)0.ida»); 
f wri  l  e  ( i  dat ,  2  jtrd.  ofl ) ; 

}  /*  for  iaamp  */ 

}/•«!» «  */ 

else  {  /*  writing  kay,  with  header  */ 
ltoa(  ldisp-  fdisp  inaroe ,  1 0 ) ; 
if  (6  >  slrlen(fname)) 

for(i=0,  cptr=header+26,  i+  +  <  6-itrten(fname);  *cptr+  +  =  ’0'); 
strcpy  (cptr.fnanw ) ; 
f  «Tite  (header.  2 .  H  E  A  D  W  O  RD  S />fl ) ; 
for(i*amp=fdi«p,  iaamp  <  ldisp ;  isamp+=nrd){ 

nrd=idatread(if3e.nrd=min(xf*iie  jdi«p-iaamp),(long)0,idat); 
fwrite(idat.2  jird,oB); 

}  /*  for  ieamp  */ 

}  /*  «u*  */ 

fc)ose(  ofl); 
relum(O); 

}  /•  wrtftle  */ 

/* 

PRMCHANGEO 

change*  GLOBALS  dynrange.  drsatirale.  colacale.  avgpow 
(pace  saving  function  for  petest 

7 

int  prmcfannge(  unsigned  ivjnt  *ag.mt  *kn){ 
int  rescale: 
rescale =0; 
switda(iv){ 

case  FDEFAULT  :  /"  restore  defaults  •/ 

•sg=l.*kn=0; 

dynrange =  DRDEFA  ULT; 

drsat  urate = D  S  D  E  FA  U  LT ; 

rescale=-l; 

break; 

case  XNOISE  : 

*kn=!*kn; 
rescales- 1; 
break; 

caae  FREQDISP  :  *sr=!*sg;  break; 
caae  DRPLUS  ; 

gamj«et(gaui_set(0X))*pow(10.0L,(MAXCOLOR*3  0L)/((MAXCOLOR-l)*20  0L))); 
colacale  *=dynr«nge/(dynrange+3); 
dynrange  +=3.0; 
break; 

caae  DRMENUS  : 

if  (dynrange  >  3.0){ 

gain-set  (gain-set  (O.D)*pow(  10. 0L,(MAXCOLOR*3  0)/((l-MAXCOLOR)*20.0L))); 

colacale  *=dynr*nge/(dynr«nge-3  0) ; 
dynrange  -=3.0; 

} 
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break; 

caae  AAPLUS  : 

gain-act  (gain-set  (0. 0)  *po  w{  10.0L,  1  0L/20.0L)); 

drsatuiate++; 

break; 

caae  AAMINUS  : 

gain-set  ( gain -aet(O.O)  *pow(  1 0.0L,- 1  .QL/20.0L)); 

draatuiate-; 

break; 

caae  FFTNCRSE  : 

xfsire  ’=  2; 
if  (xfske  >  256) 

{ 

printf("\x7");  /*  beep  to  indicate  limit  */ 
xfaiie  =  256; 

} 

break; 

caae  FFTDCRSE  ; 
xfsize  /=  2: 
if  (xfsne  <  64) 

{ 

printfC^*?");  /*  beep  to  indicate  limit  */ 
xfsite  ss  64; 

} 

break; 

}  /*  .witch  iv  */ 
return  (rescale); 

}  /*  prmchange()  */ 

/*  SCALESET 

uses  the  results  of  getscale  to  determine  scaling  for  power  spectral 
display. 

*/ 

float  scaleaet(  float  maxpow  jnt  kn){ 
maxpow  -=drsaturale, 
if(kn) 

colscale = M  A  X  CO  LO  R/ maxpo  w ; 

elae{ 

gain-act (pow{  10  0,MAXCOLOR*dynrange/(X.O*(MAXCOLOR-l ) (-maxpow /20.0)); 

colscale =(  double )( M  A  X  CO  LO  R- 1 . 0 ) /dynrangt ; 

maxpow  +  =  (double)20.0"log(  gain-set  (0.0))/lag((  double)  10.0); 

} 

ret  um(  maxpow), 

}  /’  scaleaet  */ 

r 

patest 

V 

void  patest ( void) 

{ 

int  maxval.  min  vl  ival,  is  tack .  ifl; 

ins  ktllnoiae=0,  sooag=l,  ‘sfmant ,  ’if exp  upbar=SPEKDlSP,  downbar=2; 

r 

8/6/91?  CHECK  ON  UP  AND  DOWNBAR  INDICES  AND  FREQUENCY  TRANSLATIONS 
XFSTACK  NOT  YET  IMPLEMENTED 

CHANGES  IN  XFSIZE  WILL  AFFECT  MANY  OTHER  PARAMS,  INCL.  AVGPOW 

*/ 

ins  mnvetadt[STKSZ],  mx*»tack[STKSZ].  xfatack[MAXFSIZE]; 

long  fdupUgr,  I  display,  ntansamp.  bbinaise,  iaamp.  f  curse.  1  curse; 
long  fdatack[STKSZ],  kUtadi[STKSZ],  hbinatak(STKSZ]; 
float  gainatak[STKSZ],  dacstak[STKSZ],  glob-mean,  as -noise; 
ndpage=getPage(); 
ovidmode = get  Mode)  it  i  val ) ; 

/* 

filename  (global)  filled  in  MAIN  (naig) 
or  SCRFILE: 
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•  macro  that  get*  the  filename  from  the  screm  (iplot) 
a  nufl  macro  when  compiled  as  part  of  NSIG 

V 

SCRFILE 
fcunesfdisplay=0; 
aetgmphic*(b#t  gmode ) ; 

lcurse=  ldisplay =numsamp=(fopeninp(filename  A:  ifl)- 51 2)>  >  1 ; 

writString("SIG  10/03/91  \n"  ,7,vidpage); 
xfaisesfft  -aiie(SPEKDISP); 
polar -mode(MODE-DEFLT); 
gain-aet(l.O); 

Ueekfifl.O, SEEK -SET); 

if  (HEADW0RDS<<1  !=(ival=read(ifl,header,HEADWORDS<<l))) 
errexit("  complete  header  not  found"); 
sfmant=header+122;  sfexp=header+120nampfreq=*ifmant; 
for(  ival =0  u  v»l  <  *afe  xp  aampfreq*  =  1 0  j  val + + ) ; 

*avgpow=0.0;  /*  xero  in  avgpow  aignala  first  entry  */ 
prrochange(FDEFAULT.&fdjiplay,4:killnoise); 
maxpow=getacale(ifi,  fdisplay,  ldisplay,  xfsite,  idat,  rdat, 

fc  min  val,  imaxval,  killnaiae,  hacale,  avgpow, 
Aglob-mean,  LujioIk); 
maxpo  w  sscaleae  t  ( maxpow  Jtillnoise ) ; 

/* 

now  draw  the  lines 

7 

ival=iatadi=0; 

/• 

enters  the  graphic  display  here 

starts  out  with  the  rear  graphics  cursors  selected 

left  and  right  arrows  move  the  cursor  in  sin/Je  steps  (hacale  points) 

Ctrl  key  speed  up  this  movement  5x 

home  key  selects  front  cursor,  end  selects  rec-  cursor 

pgdn  expands  such  that  selected  region  fills  the  screen 

pgup  pops  out  the  previously  selected  region 

six  levels  of  stack  are  provided 

esc  terminates  the  program 

V 

while  (EXITPROG  !=ival){ 
while)  hi  oakey  ( 1 ) ) 
bioakey(O); 

while  (lineplot(ifi.  fdisplay.  hbinsise.  minval,  maxval, 
so  nag,  Itillnoise)) 
while(bioakeyf  1 ) )  { 

if(  EXITPROG ==(  ival= bioakey  ( 0) ))  { 

writS trmigC escaped  from  screen  drawing"  ,0x4,vidp»ge); 

set  Mode  ( ov  i  dmode ) ; 

return; 

} 

else  if(prmchange(ival,ltsonag.ltkiUnoise)){ 

maxpo W  =  getsca]e(ifl,  fdisplay,  ldisplay,  xfaise,  idat, 

rdat,  liminval.  limaxval.  Itillnoise,  hacale  avgpow, 
fcglobjnean,  tem -noise); 
maxpo  w  s  seal  eae  t  ( maxpow  Xillnoise ) ; 

}  /*  if  prmrfiange..  else  */ 

}  /*  while  bioal«y(l)  */ 

/*  Make  title  */ 

golcBty(O.O); 
atrcpy(kolors  .filename); 
strcat(kolors,*  n='); 

It  oat  numaamp  Jtaiors-f  strlenf  kolors ) ,  1 0) ; 
ttrcat(kolors,"  min='); 
iloaf  min  val  Jailors  -f  strlenf  kolors ) .  10 ) ; 
strcatlkolors,"  max—") 
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itoa(  max  val  Jcolora + st  rlen(  kolors ) ,  10 ) ; 

•trcat(kolors,"  sf="); 
ltoa(sainpfre<j^ok>rs+strien  (kolors),  10); 

•treat  (kolors,”  Hz  dr="); 

iloa((inl)dynrange.kDlor»+»trl<ei(kclor«),10); 

•treat (kolors,"  aa="); 

itoa((inl)drsaturate,kolars+*trlen(hDlors),10); 

•treat  ( kolors,"  fft=" ); 
i  toa(  Kbit*  Jtolon  +*  t  rlen  (kolors ) ,  1 0 ) ; 
writStiing(kolor*,TrTLCOL,vidpage): 
if  (kiUnoise)wrilString(’  NR”,  9,  vidpage); 

/’ 

read*  the  data  between  the  cursors  (using  fdiaplay  and  hbinaize) 

•witches  to  graphic*  mode  and  clean  the  screen 

7 

fdatack(Utackj=fdisplay;  /*  save  current  location*  */ 
ldatack  pstadt  j  = lduplay ; 

if  (PSHGSTACK==(ival=getcurses(i:fdisplay.ltldisplay,  iupbar,  tdownbar, 

tfeurse,  idcurse, 

hbinaize  aonag.ldllnoise  jfl ) ) )  { 

/* 

fd  and  ldatack  (tore  current  acred  edge*  (f  and  ldiaplay  have  cursors) 
muat  aave  scaling  value*  before  incrementing  the  pointer 

7 

if  ((STKSZ-1  >  idadt)  Lit 

((fdiaplay  !=fdatack[islackj )  || 

(lduplay  !=ldatack[istack]))){ 
mnvslack{iatack]=minval; 
nu  vstack  [islack]=max  val ; 
hbinalak(iatack]= hbinaize, 
gain*  tak(iatadt]=gain-aet  (0.0); 

cl»c*tak[i»tack++]=coUcale,/*  in  ere  met*  the  »tack  pointer  */ 

} 

hbmaize=(ldispLay-fdisplay)/hscale; 

if(  (ldi^>lay-fdispla>  )%hacale) 
hbinaite++; 

lduplay = fdiaplay  4  hacale  *  hbinaize; 

maxpow=getacale(ifl,  fdiaplay.  Idisplay,  xfsize.  idat.  rdat, 
t;  min  val,  fcmaxval,  killnoiae.  hacale,  avgpow, 

L  glob  -mean.  fc«e_noiae); 

maxpow = scaleset  ( maxpo  w  JtiUnoiae ) ; 

}  /*  if  PSHGSTACK  */ 
else  if  (POPGSTACK==ival){ 

if  (0  <  iatack){ 

min  val=mnvstari[-ii  tack!; 
max  val=mx  vstack  [istack]; 
hbinai*e=hbinatak(i«tack]; 
gain-aet(  gains  tak  [is  tack] ) ; 
colacate=dacatak[iatackj ; 

} 

fdiaplay =fdatadt[iatack] ; 

ldiaplay = ldatack  [iat  ack  ] ; 

}  r  if  POPGSTACK  •/ 

elae{ 

if  (WRTFILE==ival) 

if  ( wrtfilef  ifl  JduplayJdi*play)){ 
gotcay(O.O); 

writStrkigC  error  writing  Rle’Dxf.vidpage); 

} 

/* 

geteuraea  returns  the  cursor  positions  in  fdiaplay  and  ldisplay 
muat  recower  true  display  positions  from  the  stack  if  no  change 
desired,  as  in  WRTFILE  and  FRQD1SP 

7 
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fdisplay =fdstack[istack] ; 
ldisplay— ld»lack[i»tack]; 
if  ( prmchange(  i  val  ,&»onag  ikillnoiae) )  { 

get*caie(ifl,  fdiiplay,  lditplay,  xfiiie,  idat,  rdat, 

i:  nun  val,  i:  mar  val,  killnoiae,  hacaie,  avgpow, 
ti glob-mean  &u_noi ae); 
marpow  =acaleset  ( maxpo  w  Jrillnoiae ) ; 

} 

} 

if  (numump  <  (ldispiay=fdisplay  +  (long)hacale*hbinai»e)) 

1  display = n  urnaamp ; 

}  r  while  EXITPROG  1 — i val  •/ 
cloaefifl); 

aetMode(ovidmode); 

return; 

}  /*  pateat  •/ 
main(argc,argv) 
im  argc;  char  “argvQ; 

{ 

if  (1  ==  argc) 

errexit("uaage:  naig  <filename>  dr  aa  [e|v]“); 

•trcpy  ( filename  .argv  [l  ] ) ; 

if  (V  ==  (argv[4](0]  |  0*20)){  /*  binary  or  convert!  to  lower  caae  */ 
vacate  =  EVSCALE; 
batgmode  =  EBSTGMODE; 
botline  =  EBOTLINE; 

}  /*  default,  are  VVSCALE,  VBSTGMODE.  VBOTLINE  */ 
vacale-  =  16; 

if  (!(dynrange  =  atof(argv[2]))) 
dynrange  =  DRDEFAULT; 
if  (!(dreaturate  s  atoi(argv[3)))) 
draaturate=DSDEFAULT ; 
pete*t(); 
retum(O), 


} 
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dac.c 

This  function  reads  and  tend*  the  kay  file  data  to  a  Canetics  board  model 
PC-DMA.  The  D/A  capabilities  of  the  Canetic*  board  are  used  to  'pUy'  the 
kay  file  while  the  time  atria  and  time-frequency  distribution  are  displayed 
by  SIG. EXE. 

The  following  functions  were  taken  from  the  Model  PC-DMA  User  Manual, 

Canetics,  Inc. 

memreqO 

init  _dma_ad( ) 

init-dma-da( ) 

init-d  .a.da-cont  ( ) 

init-board() 

kilLtimer() 

init  -timer  da() 

The  following  function  was  taken  from  Graphics  Programming  in  C,  Roger  T. 

Stevens,  M&T  Publishing,  Inc.,  Redwood  City,  California,  1988. 

plot() 

7 

#include  <Libl80h> 

#indude  <stdio  h> 

#include  <stdlib  h> 

#include  <dos.h> 

#include  < alloc  h> 

#include  <conio.h> 

#include  <fcntl.h> 

#include  <stat.h> 

# include  < string. h> 

#include  <io.h> 

#include  <fcnil  h> 

♦  include  <math.h> 

♦define  CONTROLA  /*  add+lM  */  15876 
♦define  CONTROLB  /*  add+2’4  */  15880 
♦define  OFFSET-ADD  /*  add+lCT4  7  15912 
♦define  RESET -ADD  /*  add+ll*4  */  15916 
♦define  DMAMODE-ADD  11 

♦define  DMA-M ASK-ADD  10 

♦define  CLRPTR-ADD  12 

♦define  AD-TIMER-MODE  122  /*  ad  is  timer  ^1  in  this  mode  */ 

♦  define  DA.TIMER-MODE  52  /*  and  da  is  timer  ^0  "/ 

♦define  TIMER-AD  15925 

♦  define  TIMERJDA  15924 

♦  define  TIMER-ADD  15924 

♦  define  TIMER-MODE-ADD  15927 

♦define  STA  15904 

♦define  STD  15908 

♦define  TOGGLE-DAC-ADDRESS  15888 

♦define  UH-SEG(p)  ((  FP-SEG(  p  )  4-  ((FP.OFF(  p  )  )  >>  4)  )  L  OxFOOO) 

♦define  UH-OFF(p) 

((int)((((FP-SEG(p)-t-((FP.OFF(p))>  >4))l:OxOFFF)<<4)+(FP  -OFF(p)fcOxOOOF))) 

typedef  struct  settings  *pseU; 
struct  setting  setsl; 
pacts  sets; 
char  curfibv*nie[30j; 
unsigned  long  data-sjxe; 
ini  mem  -allc=0; 

ini  far  *memreq(void)  /*  Allocate  two  pages  of  memory  and  choose  */ 

{  /*  one  complete  page  */ 
ini  tempseg; 

void  far  *death-lo-nonbeliever*; 

death -to-nonbelie  vers  =  farmed  ocl  0x2  OOOOL); 

if  ( death -to-nonbe  Lie  vers  ==  NULL) 

{ 

printf("Corelrft  =  Hid  Need  %ld\n\r\  corelrftl).  Cb<20000L  ); 
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putchar(l); 

} 

tempseg  =  L'H-SEG(death-lo-nonbelievers); 
if  (UH-OFF(death-lo-nonbelievax)  !=  0) 
tempseg  +=  0x1000; 
return  (int  far  *)MK-FP(tempseg,  0); 

} 

void  Lai t^imA-Ad( unsigned  couni ,  unsigned  *eg) 

{ 

outp(DMA-MASK_ADD,  0x05); 
outp(DMAMODE-ADD,  0x45); 
outp(CLRPTR_ADD,  0); 
outp(2,  0); 
outp(2,  0); 

outp(3,  count  L  OxOOFF); 
oulp(3,  count  >>  8); 
outp(130,  teg); 
outp(131 ,  seg); 
outp(DMA-MASK_ADD,  1); 

} 

void  init-dma-da(  unsigned  count,  unsigned  seg) 

{ 

outp(DM A -MASK-ADD,  0x07); 
outp(DMAMODE_ADD,  0x4B); 
outpfCLRF^R -ADD,  0); 
outp(6,  0); 
outp(6,  0); 

outp(7,  count  &  OxOOFF); 
outp(7,  count  >>  256); 
outp(130,  seg); 
outp(131 .  seg); 
outp(DMA-MASK-ADD.  3); 

} 

void  init-dma  da  root  (unsigned  seg) 

{ 

outp(DMA-MASK_ADD,  0x07); 
outp(DMAMODE_ADD,  0x5B); 
outp(CLRPTR_ADD,  0); 
outp(6,  0); 
outp(6,  0); 

outp(7,  OxFF); 
outp(7,  OxFF); 
outp(130,  seg); 
outp(131.  seg); 

outp(DMA -MASK -ADD,  0x03), 

} 

void  uut-board( struct  settings  ‘sets) 

{ 

outp(RESET_ADD,  0); 
outp(OFFSET_ADD,  seta->offset); 

outpfCONTROLA,2O0);  /*  PC-DMA  Control  register  A:  11010001  */ 
outp(CONTROLB,154);  /*  PC-DMA  Control  register  B:  10011010  ’/ 
outp(RESET_ADD.  0); 

} 

void  kill -timer)  struct  settings  ’sets) 

{ 

outp(R£SET_ADD,  0); 

outp(CONTROLB,  15  +  16*sets->  gambits  +  128); 
outp(RESET-ADD,  0); 
if  (sets- >dac2fiag  =  —  2) 

outp(TOGGLE  JDAC-ADDRESS,  0); 
outp(TIMER-MODE-ADD,  SO); 
outp(TIMER-MODE-ADD,  114); 
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outp(TIMER-MODE-ADD,  178); 
injt.dma-da(0,  sets->seg); 
init-dma-ad(Q.  sets->seg); 

} 

void  init-limer-da(unaigned  spaces) 

{ 

outp(TIMER_MODE-ADD,  DA-TIMER-MODE); 
outp(TLMER-D A,  tpAces  L  OxOOFF); 
outp(TLMER-DA,  spaces  >>  8  ); 

} 

void  plot  (ini  x.  ini  y,  int  color) 

{ 

unsigned  int  offset; 

ini  dummy,  mask; 

char  far  *roem-addreas; 

offset  =  (long)y*80L  +  ((long)x/8L); 

mask  =  0x80  >>  (x%8); 

-ES  =  Ox  A 000, 

_BX  =  offset; 

.CX  =  color; 

-AX  =  mask. 

asm  MOV  AH.AL 

asm  MOV  AL.08 

asm  MOV  DX.03CEH 

asm  OUT  DX.AX 

asm  MOV  AX.0FF02H 

asm  MOV  DL.0C4H 

asm  OUT  DX.AX 

asm  OR  ES;[BX],CH 

asm  MOV  BYTE  PTR  E5:[BX],00 

asm  MOV  AH. CL 

asm  OUT  DX.AX 

asm  MOV  BYTE  PTR  ES  [BXj.OFFH 

asm  MOV  AH.0FFH 

asm  OUT  DX.AX 

asm  MOV  DL.0CEH 

asm  MOV  AX.0003 

asm  OUT  DX.AX 

asm  MOV  AX.0FFO8H 

asm  OUT  DX.AX 

} 

void  wtfile.cont(psets  sets,  mt  curftl) 

{ 

unsigned  ii,  jj,  tempeount.  ckpoint=0,  last,  iark,  buffcr(0x4000],  loc=0; 

unsigned  far  *blodi,  far  "uptr, 

unsigned  near  *nptr; 

unsigned  long  chunk,  mAx.pt,  point, 

double  cununc; 

block  =  MK -FP(set*  >seg  <<  12,  0).  /*  Pointer  to  allocated  memory  9 / 

kiU-limer(seu); 

max  _pt  =  daia-site: 

cursanc  =  (double  )640/(doublf)data  si*e; 
plot(0, 215,1  5); 

if  (d  At  a  site  <  0x10000)  /*  If  <  one  page,  read  in  all  dAts  and  9  f 
{  /•  then  cycle  to  PC-DMA  9 / 
ploMO, 215,0). 

point  =  ceil ( cun unc*( double )(mAx -pt-dsts  >ite)); 
plot(  point, 21  5,1 5); 
uptr  =  block; 

»hile(dat  A-Aise  >  0x4000) 

{ 

reAd(curfU.  buffer,  0x4000); 
data-size  -=  0x4000, 
loc  +  =  0x4000; 
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ckpoint  +=  64; 
nptr  s  buffer; 

while  (FP-OFF(uptr)  <  loc) 

*uptr+  +  =  *nptr++  <<=  4; 

} 

plot(  point, 215,0); 

point  =  ceU(cur»_inc*(double)(max-pt-dat*  »ite)); 
plot) point, 21 5,15); 
re*d)curfii,  buffer,  d*l*  lire); 
loc  +  =  d*t»_*ite; 
dcpoint  +=  64; 
if  (ckpoint  r- =  319) 
ckpoint  =  63; 
nptr  =  buffer, 

while  (FP-OFF(uptr)  <  loc) 

*uplr+  +  =  *nptr+  +  <<=  4; 

•wme  =  *uptr 

while  (FP.OFF(uptr)  <  loc  -  dtU  tite  +  0x4000) 

*uptr+  +  =  unx; 
ini  t -board)  MU); 

init  dm  A  d»  <-nnl(»»U.  >teg); 
init-timer-d*(»et*->wait»); 

outp(CLRPTR_ADD,  0); 
inp(6); 

tempcount  =  inp(6); 
while  (lempcour*  <  ckpoint) 

{ 

inp(6); 

tempcount  =  inp(6); 

} 

} 

if  (daia-*ize  >  0x10000)  /*  If  >  one  page,  fill  one  page  and  then  */ 
{  /*  fill  quarter  page  u  DMA  paaae*  proceed*  */ 
uptr  =  block: 
for  (ii=0 ii<4ji  +  +  ) 

{ 

read(curfil.  buffer,  0x4000); 
datafile  -=  0x4000; 
nptr  =  buffer; 
loc  +=  0x4000; 
if  (loc  ==  0) 

loc  —  Oxfffe; 

while  (FP-OFF(uptr)  <  loc) 

*uplr++  =  *nplr+  +  <<=  4; 

"uptr  =  “nptr  <<=  4; 

} 

uptr-f  + ; 

init_board(*eta); 

ini  t.  dm*  riaront(»eU-  >*eg); 

init_un»er.d*(*et*->wait«); 

outp)CLRPTR_ADD,  0); 

plot)0, 215.0); 

point  =  ceil  (cur* -inc*(double )( max -pt-data.*ite-Qx  10000)); 

plot) point.  21 5.1 5), 

inp(6); 

tempcount  =  inp(6); 
while  (tempcour*  <  63) 

< 

inp(6); 

tempcount  =  inp(6); 

} 

for  (ii  =  l;ii<*et*->num-of-pag;ii  +  +  ) 

< 


clrpoint  =  127; 
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SIG  —  D  ACC 


chink  =  0x4000; 
for  (jj=l;  jj<=4;  jj++) 

{ 

read(curfil,  buffer,  0x4000); 
data  »ite  -=  0x4000; 
plot(poinl, 215,0); 

point  =  ceil(ct«x  _inc*(double)  (max -pt-data-site- 0x1 0000)); 
plotfpoim, 215,15); 

□ptr  =  buffer; 

while  (FP.OFF(uptr)  <  chunk) 

*uptr+  +  =  *nptr++  <<=  4; 

*uptr  =*nptr  <<=  4;  /*  ch  */ 
inp<6); 

tempcount  =  inp(6); 

while  {  tempcount  <  ckpoint  ) 

{ 

tnp(6); 

tempcount  =  inp(6); 

} 

clmmk  +=  0x4000; 
if  (chunk  ==  0x10000) 
chunk  -=  2; 
ckpoint  +=64; 
if  (ckpoint  ==  319) 
ckpoint  =  63; 

} 

uptr++; 

} 

last  =  ceil((double)data_aiie/0x4000); 

ckpoint  =  127; 

chunk  =  0x4000; 

for  (jj=  1 ;  jj<last;  jj  +  +  ) 

{ 

read(curfil.  buffer,  0x4000); 
data_siie  -=  0x4000, 
plot(poinl,  215,0); 

point  =  ceil(curvunc"  (double)  ( max  _pt-dala-sixe-0x  10000)); 

plot(point.215.15); 

nptr  =  buffo-; 

while  (FP-OFF(uptr)  <  chunk) 

*uptr++  =  *nptr+  +  <<=  4; 

*uptr  =  *nptr  <<=  4; 
inp(6); 

tempcount  =  inp(6); 

while  (  tempcount  <  ckpoint  ) 

{ 

inp(6); 

tempcount  =  inp(6); 

} 

chunk  +  =  0x4000; 
if  (chunk  ==  0x10000) 

chunk  .=  2; 
ckpoint  +  =  64; 
if  (ckpoint  ==  319) 
ckpoint  =  63; 

} 

plot!  point, 2 15.0); 

point  =  ceil(cur»  -inc*(double)(max-pt-0xl0000)); 

plot(point.215,15); 

readicurfU.  buffo-,  data  sire); 

nptr  =  buffer, 

while  (FP.OFF(uptr)  <  chunk  -  0x4000  +•  data  site) 

*uptr++  =  *nptr++  <<=  4; 

*uptr++  =  *nptr-++  <<=  4; 


SIG  —  DAC.C 

inp(6); 

tempcount  =  inp(6); 

while  (  tempcount  <  drpoint  )  /*  Watch  DMA  and  wait  until  data  */ 
{  /*  haj  all  been  tent  to  board  */ 
inp(6); 

tempcount  =  inp(6); 

} 

plot(  point, 215,0); 

point  =  ceil(cur».inc*(ck>uble)(max-pt-0x4000)); 
plot(point,215,15); 
dt  point  +=  64; 
i f  (depoint  ==  319) 
ckpoint  =  63; 
inp(6); 

tempcount  =  inp(6); 
while  (  tempcount  <  255  ) 

{ 

inp(6); 

tempeouni  =  inp(6), 

} 

plot(  point, 21 5.0), 

point  =  ceil(curt-inc*(<iouble)(mAX_pt-0x8000)); 
plol(point.215.15); 
if  (ckpoint  ==  63) 
ckpoint  =  255; 

eUe 

ckpoint  -  =  63; 
inp(6); 

tempcount  =  inp(6); 

while  (  tempcount  <  ckpoint  ) 

{ 

inp(6); 

tempcount  =  inp(6); 

} 

plol( point. 215,0); 

point  =  ceil(cun-inc*(double)(max_pt-0xa000)); 
plot(  point.  215, 15); 
if  (ckpoint  ss  63) 
ckpoint  =  255; 

el*e 


ckpoint  -=  63; 
inp(6); 

tempcount  =  in  p  ( 6 ) ; 

while  (  tempcount  <  ckpoint  ) 

{ 

inp(6); 

tempcount  =  inp(6), 

} 

> 

lull-timerfeeti); 
plot  (point  .215,0); 

} 

vad  >tart-dac(int  ifile,  long  'begin,  long  *end.  long  aamplej-ate. 
Boat  teal) 

{ 

int  far  *M  J’tr; 
int  *tmant.*texp; 
char  he*d»(512]; 

outp<RESET-ADD.Oxda);  /*  Detect  Canetict  board  */ 
if  (inp(  RESET -ADD)  -  =  0)  /*  0  indicatet  board  it  preaent  */ 

{ 

tett  =  ittettl:  /*  Get  the  t  true  lure  tett'  variable*  */ 
if  ('mem-allc)  /*  Allocate  memory  only  rnce  */ 

{ 
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M-Ptr  =  memreqO; 

*ets->seg  ss  FP-SEG(  M-Ptr  )  >>  12; 
mem-aUc=l; 

} 

if  ('begin  ==  0)  /*  Move  to  start  in  kay  file  '/ 

lseek(ifile,512,SEEK-SET); 

else 

l*eek(ifile,'begin'2  +  512.S£EK-SET); 
data-size  =  2'('end  -  'begin); 
sample-rate  *=  seal; 

sets- >  waits  =  floor(le6/sample-rate  +  0.5); 

sets-  >oum_of_pag  =  floor)  (double  )data-#ise/0»  10000); 

sets-  >  daiasne  =  1; 

sets-  >  miine  =  0; 

sets-  >  runs  top  =  0; 

sets- >gambsts  =  1; 

sets- >  inmode  =  2; 

sets-  > offset  =  128; 

sets-  >  Alternate  =  1; 

sets-  >  uBngx  mux  =  0; 

sets-  >dac2flag  —  1; 

injt -board  (sets); 

kiU-timer(sets); 

wtfiie-cons(sets.tfile); 


SIG  —  DRAWHLS.C 
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/*  byte  mash  always  Orff,  color  in  di,  mask  in  si,  j libel  replaces  noxor  */ 
/*  bx:  memory  location;  dx,  ax:  genval  purpose  register*  */ 

/* 

DRAWhlN 

ASSUMES  SMALL  MEMORY  MODEL 

draws  a  single  pixel  votical  line,  located  at  homontal  position  x 
vertical  bounds  of  the  line  are  yl  and  y2 

xor  function  can  be  selected  with  colcr=  color*  16,  overwrite  otherwise 
c change  s  increment  for  color  pointer 
*/ 

# pragma  inline 

void  drawHLn(inl  y,  int  xl,  int  x2.  char  color){ 

/* 

exchange  yl  and  y2  if  yl  larger  than  y2 

V 

if  (xl  >  x2){ 

asm  tnov  bx.x2 
asm  mo v  ax, il 
aam  mov  x2.ax 
asm  mov  x  1  ,bx 
} 

/• 

any  use  of  egacgx  requires  bracketing  by  the  six  port 
instructions  bracketing  the  for  loop  below 
see  Dwyer  el  al  micro/systems  July 88  20-34 

*/ 

/*  setting  up  ega  for  write  mode  2  */ 
outportb(Ox3ce.5); 
ou  t  portb  ( Ox  3cf .  2 ) . 

/*  set  cl  =  xl  mod  8  */ 
asm  mov  ax.xl 
aam  mov  cl.al 

/*  set  ch  s  xl  div  8  ranges  from  0-80  */ 

aam  aar  ax.l 
asm  sar  ax.l 
aam  aar  ax.l 
asm  mov  ch.al 
/’  bx  =  80  *  yl  */ 
aam  mov  ax.y 
aam  sal  ax.l 
asm  sal  ax .  1 
aam  sal  ax.l 
aam  sal  ax.l 
aam  mov  bx.ax 
aam  sal  ax.l 
aam  sal  ax.l 
asm  add  bx.ax 
/*  bx  =  80  •  y  +  x/8  */ 
asm  xor  ax  .ax 
mm  mov  al.dh 
aam  add  bx.ax 
/*  bx  *  addrea  offset  */ 

/'as  ega  ram  startmg  address  */ 
aam  mov  ax.OxaOOO 
aam  mov  a  ax 
/*  set  di  =  color  */ 
aam  mov  di.  color 
aam  and  di.Oxlf; 

/*  setup  complete,  remaining  item#  vary  */ 

/*  now  plot  first  partial  byte,  seromg  the  initial  bits,  put  mask  in  ai*/ 
aam  mov  al.Oxff 
aam  and  cl.  0x7 
aam  js  takip 
aam  shl  ai  d 
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SIG 


DRAWHLS  C 


asm  ihr  »J,d 
itkip : 

aam  xor  ah. ah 
aam  mov  u.ax 

/* 

select  bit  maak  register 
output  bit  maak 
back  to  select  register 
function  select  register 
reset  bits  3  and  4 
check  for  xor  function 
set  bits  3  and  4 

7 

asm  mov  Ax  0x3ce 
aam  mov  el, 0x8 
asm  out  dx.al 
«im  me  dx 
aam  mov  ax, si 
aam  out  dx.al 
aam  dec  dx 
aam  mov  al.0x3 
aam  out  dx.al 
aam  inc  dx 
aam  xor  al.al 
aam  test  dj.OxlO 
aam  jx  injlp 
aam  mov  al.OxlS 
initp: 

aam  out  dx.al 

/*  maak  and  function  selected  */ 
aam  mov  ah.es  [bx| 
asm  mov  ax  .di 
aam  mov  es  fbxj.al 

/*  now  prepare  for  all  whole  byte  plots  */ 

/*  compute  x2  div  8.  and  subtract  xl  div  8  */ 
aam  mov  ax.  x  2 

aam  sar  ax.l 
aam  sar  ax.l 
asm  sar  ax.l 

aam  tub  al.ch 
aam  dec  ax 

asm  jx  Leal 

/* 

save  loop  count  in  cx 
select  bit  maak  register 
output  bit  mask 

7 

««ni  mov  cx.ax 
aam  mov  dx.0x3ce 
— m  mov  ai  0x8 
a am  out  dx.al 
asm  inc  dx 
aam  mov  ex.Oxff 
aam  out  dx  .al 
/*  loop  •/ 
ploop: 

aam  inc  bx 

aam  mov  ah  .ea  fbx] 

aam  mov  ax.di 

mov  es:[bxj.al 
asm  loop  ploop 

/*  now  prepare  last  partial  byte  */ 

/*  set  d  =  x2  mod  8  */ 

last 


SIG  —  DRAWHLN.C 


Mm  mov  ax.xj 

Mm  mov  cl. 0x7 
Mm  ud  d.al 
Mm  me  bx 
Mm  mov  «J.0xfT 
Mm  xor  cl, 0x7 
Mm  jx  L»kjp 
Mm  <hr  ad  d 
Mm  (hi  al.d 
Ukip: 

Mm  mov  >1  ax 

/* 

(elect  bit  mMk  register 
output  bit  mMk 
back  to  select  register 
function  (elect  remitter 
met  bite  3  and  4 
check  for  xor  function 
(et  bita  3  and  4 

7 

asm  mov  dx.OxOce 
Mm  mov  al.0x8 
aim  out  dx.al 
asm  inc  dx 
Mm  mov  ax.ai 
asm  out  dx.al 
Mm  mov  ah  et  [bx] 
aim  mov  ax  di 
asm  mov  e«  [bx],al 
outportb(0x3oe.5); 
o  u  t  portb  ( Ox  3<d ,  0 ) . 
ou  t  portbl  Ox  3ae .  8 1 ; 
outportb(Ox3cf.OxfD. 

/*  returning  ega  to  write  mode  0  •/ 
}  /*  draw v In  */ 
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SIG  —  FFTFUNCS.C 


#  define  MAXBOOST  -15.0 
*-  'Jude  "aagbla-h* 

r 

thii  file  con  air*  all  of  the  A  rallied  code,  «n«  functiraia 

that  call  the  fft  routatn 

it  alao  contain  a  few  Matic  global  variable!  that  are  particular  to 
the  fft  and  related  lubrotaiaae 

V 


/*  MUST  CHANGE  THE  MAXBOOST  CONSTANT  AND  ITS  PURPOSE  •/ 
extern  void  far  fwil 

rfft  (float  far  ‘dat,  iat  far  ‘exp,  iat  far  snorm1  far  *»d); 
extent  void  far  fieril 

polar<Baat  far  *idata,  float  far  *odal,  loot  far  *adataaa,  iat  far  ‘anode); 


r 


hainm( float  far  ‘rdata.  Boat  tar  ‘idat,  Lr> 


GLOBALS  uaad  by  Micro  Way  87  routine 

V 


etatic  iat  nortnaliae,  poerj;  /•  uaad  by  rfft  */ 

•tatic  iat  node,/*  6  foe  naig,  ipkot;  1  for  kayuau  •/ 
italic  long  dataiae;  /*  uaad  by  polar  •/ 

•tatic  float  gala;  /•  uaad  by  rfft  •/ 


/’ 

theae  function  iaitialke  aad/or  report  GLOBAL  variabln 
uaad  iB  Micro  Way  routines 

GLOBALS  MUST  NOT  BE  MODIFIED  ANYWHERE  ELSE  IN  FFTFUNCS,  AND 
THESE  PUBLIC  FUNCTIONS  ARE  THE  ONLY  ACCESS  FOR  OTHER  SOURCE  FILES 


iat  ffLaite(iot  xf){ 

aormalac  *  2;  /*  rfft  •/ 
if(xf){ 

datain  *  (xf>>l);  /•  polar  •/ 
po«3  ■  0;  /*  rfft  */ 
while(  xf>  > + + (poo  2 ) ); 

(powJH 

} 

return!  ( iat  )datai«e  <  <  1 ) ; 

>  r  Ajim  ’/ 

iat  polar -mode  (iot  pja){ 
if(p^o) 


retuxa( mode); 

> 

Boat  gaia-aet(  double  ga){ 
tf(ga) 

gain*  ( float  )pt; 
ratumf  gaia); 

r 

FLDEMEAN 


cMpeta  aad  nibtraru  out  the 
at*  FOR  loop  require!  poiter 

V 


mean  at  a  float  array 

arithmetic  aad  logical 


raid  Bdamaaa(float  *fdata,  iat  iaiae){ 

Beat  “fptr;  double  mnn 

focfflxrmfdat nfiaiae.  nmaaoO.O;  -fptr>*fdata;  maea+**fptr); 

forffrxrolditu  iaiae;  -fptr>*fdata,  ‘fpa-.snaaa); 

}  r  Bdeaaaaa  •/ 

r 

IMEAN3FP 


detnaaaa  integer  data  in  idat,  placet  it  la  rdat 

V 

void  imeaajfpfml  ‘idat.  float  *rdet,  iat  iaiae){ 


SIG  —  FFTFUNCS.C 
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int  "iptr, 
loaf  imean; 
float  "fptn 
double  dmeu; 

for( nnean=0.iptr=]cUt-fitiie;  -ipcr>=idat;  unean+  =  "iplr): 
dmean  =  (double)iroean/(doufele)i*jK; 

for((ptr=rdat  + iaise.iptr=idjt-f mu;  fptr>rdat;  *-fplr=*-iptr-dmean); 

} 

/• 

MAKEpSPEC 

calls  canard  liamming  window,  canned  fft  routme  for  real  data 
converts  complex  results  to  maf ,  phase  and  returns  in  rdat 
GLOBALS  INITIALIZED  IN  INITFFT() 

7 

void  mak*Pspec(float  *rdst,  int  isisc){ 
hamm(rdat,  rdat,  Xcpowl); 

*(rdat+ (isi«e>>l))  +*  0.001;  /"  prevents  underflow  in  log  power  */ 
rfft(rdat,  k pow2,  normalise.  fcgain); 
polar  (rdat.  rdat.  fedataise,  Xtmode); 

}  /"  makePspec()  •/ 

/• 

NORMpSPEC 

subtracts  the  power  spectrum  values  in  svgpQ  from  rdatfl 
note  that  rdat  has  PS  values  in  even  entries  only 
avgp  has  PS  values  in  every  entry 

V 

void  noimP spec ( float  'rdat,  float  "avgp,  ins  iaisc){ 
float  *fsif; 

for(fcigzrdat-f  iaise,  avgp-f  s(isise>>l); 
fa.  f>  rdat,  "(faig-  =  2)-**-avgp); 

} 

/* 

COLRpSPEC 

converts  power  spectrum  values  in  rdat  to  integer  color  values  in  kolrs 
relying  on  pointer  arithmatic  and  logical  comparisons  for  iptr.  kolrs 
this  version  allows  fft  sixes  other  than  256 

V 

void  colrPspecffloat  "rdat,  char  "kolrs,  int  iaise,  float  colscale){ 
char  c; 

ins  "iptr,  ncelb,  k.  j; 
float  "fptr 

nceUsnSPEKDISP/iaiae; 

/"  iptr  is  placed  SPEKDISP  bytes  beyond  the  beginning  of  kolrs  "/ 
fori  fptr  -l  rdat  -f  isise ,  iptr=  kolrs -r  SPEKDISP;  iptr  >  kolrs;  ){ 
k  *  (ins)(colncale"  *(fpti**2)); 
k  s  min(MAXCOLORjnax(OJt)); 
farijwO;  j++<ncells;  "-iptr=(k|(k<<8))); 

> 

}  /•  COhPspWcO  "/ 

r 

SDIFPCMP 

V 

ins  adiifacnpf coest  void  "reel,  const  void  "rec2){ 
float  fl,  *f2; 
fl  *  reel-.  f2  *  reel, 
if("fl>*f2) 
return(l); 
elaa  if  ("fl<*f3) 
return(-l); 

elaa 

retvn(O); 

} 

/* 

GETSCALE 
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SIG  —  FFTFUNCSC 


7 

Boat  gettCAle(int  ifp,  loaf  fdirp,  loaf  Iditp,  ml  i/vu, 

mt  ideiQ,  Boot  rd*t[],  ini  *imin,  int  *i mex, 
in t  momt ,  int  htcele,  Bom  *vgpQ, 
floet  float  *t«moLM){ 

typedef  itruct  {float  M,  int  mm;}  etblk; 
etblk  blk[HSCALE]; 

int  *iptr,  iold,  i,  diff.  jmm,  jmin,  n  block*,  nttop,  imtna; 
int  firethiok.  Lattbbk; 

floet  lunaq,  otui,  *fptr,  *dptr,  mrpow,  global-mean; 
double  xsJDeu; 
lon(  Ukip.  fplace; 
gaineet(l.O); 

memeet(blk.0.aii*ctf(eeblk)*HSCALE); 

laeek(  ifp, (fdiep-t  HEAD  WORDS)<<l, SEEK-SET); 

“imax  s  -(*imin  s  MAX  INT); 
for ( global-me an  —  mtut  *  imaxat  =  i  =  0, 
nblockeaminf  hteeie ,  (Iditp-  fdUp )  /  xfti  te ) . 

Ukip  *  max(0, (Iditp  •  fdi*p)/hacal*-xfai*e); 

i  <  n block*,  !++){ 

fpleceatell(ife); 

nump  =  id*Qtad(  Jp.xfaixeJaiupjdat); 

/*  file  handler.  #itemt.  byte*  to  akip  After  reed,  *deatination,  ret  #inU  */ 

for(rum*q=0.0. global-mean -f  sioldsjmAxsjminz^iptraidai+naainp-l); 
-iplr>  =  idat,  ){ 
globe!  mean+ =*iptr: 
dilT  a  ’iplr-iold; 
lumq  +*  difPdUf; 
if("ipcr  >  jmu) 
jmax  a  *iptr, 
else  if  (*iptr  <  jmin) 
jmin  a  *iptr, 
iolda°iptr. 

} 

if  (“imax  <  jmex) 


if(*lmin  >  jmin) 

"imiw  S  jmin; 

blk[blk[i]  jxun  *  i]e»  *  tuzneq; 

if(maam<tumeq){  /*  needed  for  iwekng  color*  during  toom  */ 
mtinsnuaq;  /*  when  full  qaort  it  not  performed  */+ 


}  /•  for  i  0  to  nbiock»-l  •/ 

*  g-meens  globeLmeen /(( double  )nblockt  *  xiute ) ; 

Ukip  +■  xfaiae;  /*  Ukip  beromre  reed  inlervel  */ 

/*  poler  jnode(*treight  Amplitude)  */ 

if  (*avgpaa0.0){  /•  only  done  once  for  the  Ale  */ 
qeort(blk,rf)lr»rlre  eieenf(ed)lk)  ediffanp); 
memeet(*rgpjO,(l+(xf*i**>>l))*ei*eof(Boet)); 

firethlok  a  nblockt/10; 
if('(lmt  blckan  block*/ J)) 

Ltetblak-f-f  ;/*  10-30%  */ 

forfiaflntblok;  i<lettbldt;  /•  i++  in  loop  body  */){ 
em  meen-f  »blkfi)  **• 

Ueek(i^,(HEADWORDS+blk(i++].ni*n*Ukip)<<l  SEEK-SET); 
neemp  *  idetreedf ifp jrfuee,( long )0*det); 
im*en2fp(  idet  jdet ,  xfeUe ) ; 
mekePcpecjrdet  .xfttte); 

forfdptraAvgp-f  xfaiee/3,  fptr  a  rdal-f  xfeiee;  (fpa^  =  2)>*Tdai; 
*-dptr  -fa  *fptr); 


SIG  —  FFTFUNCS.C 


/* - compute  for  pow«  ipectra - -*/ 

/• 

THIS  NEEDS  WORK:  BETTER  MAXBOOST  VALUE 
if  polar -mode  react,  need  to  compute  201n(*dptr/iisneq) 

*/ 

for(dptr=avgp-frfaiae/2,  lumaq  s  Laatblok-finthJok: 

-dptr>=avgp;  ) 

if(MAXBOOST  >  (*dptr  /=  lunaq)) 

•dptr=MAXBOOST; 

}  /•  if  (*a»B>s=0,0)  */ 

/* - goto  maximum  block - */ 

/*  polar  .mode)  20ln( amplitude))  */ 

if  ((long )-l a slaeek(ifp,(lom)(fdi»p-i-imar«a*Ukip-l-HEADWORDS)<<l, SEEK-SET)) 
errexit  ("error  in  aeekiag  maximum  •umaquerea  block"); 
i  dal  mad)  ifp,  ( int )  xfaiae ,( loog)0  jdat ) ; 
imenn2fp(idat.rdat,xfii»e); 
mahcP«pec(rdal,rfai»c); 
if  (xnoiae) 

nortnP«pec(rdal,avjpjcf»ixe); 
for(mrpo»=*(fptr=rdal).  isl; 
i++  <  (xf*ise>>l);){ 
if(*(fptr+=2)  >  mxpo*) 
mjtpo  w  s  *fptr, 

} 

return)  mxpo*); 

}  /•  ffcale  */ 
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SIG  —  LIB180.H 


(truct  Kttings 
{  unsigned  mime; 
unsigned  runsteg , 
unsigned  gambits, 
unsigned  Alternate; 
unsigned  dscssixe: 
unsigned  inmode; 
unsigned  (eg; 
unsigned  offset; 
unsigned  waits; 

unsigned  long  num-of-con; 
unsigned  dacSflsg; 
unsigned  nmn-of-pag, 
iimifiw d  uflinfxxnax; 

}! 

void  Init  »d(nn«y..d  count,  unsigned  Kg); 
void  init.dma-da( unsigned  count,  unsigned  Kg); 
void  ini  t.dma.  ad -coal  (unsigned  seg); 
void  ini  t-dma-da_coat(  unsifted  Kg); 
void  tnii-board(  struct  ki tings  *Kts); 
void  killaimer(stnsct  settings  *Kts); 
void  injtaimer-de( unsigned  spaces); 
void  inn  timer  sd( unsigned  spaces); 

void  «Titesingle( struct  Kttmp  *Kls,  unsigned  msb,  unsigned  Isb); 

void  readfile<  struct  settings  *ku); 

void  wntefile( struct  settings  *Kls); 

void  readfile-cotn  ( struct  Kttings  *Kts.  char  'filename); 

void  wntefiJemonl (struct  Kttings  *Kts.  chsr  'filename); 

ml  status! mt  print); 

unsigned  readsingle( struct  Kttings  *Kts); 

# define  G.VERSION  '1.VT 


SIG  —  NDRAWVLN.C 


r 

DRAW  VI N  11/24/88 

ASSUMES  SMALL  MEMORY  MODEL 

draw*  a  tingle  pixel  vvtical  line,  located  at  horoontal  poaition  x 
vertical  bound*  of  the  line  are  yl  and  y2 

xor  function  can  be  aelectcd  with  color=color  + 16,  overwrite  otherwiae 
!thia  function  ii  presumed  to  be  aelected  for  the  entire  Line! 
cchange  =  increment  for  color  pointer 

*/ 

#  pragma  inline 

void  drawVLnfint  x,  ini  yl,  ini  y2,  char  *colcr,  int  cchange ) { 

/• 

exchange  yl  and  y2  if  yl  largs  than  y2 
*/ 

if  (yi  >  y2){ 

a am  mov  bx,y2 
a am  mov  tx.yl 
aam  mov  y  2  ax 
aam  mov  y  1  ,faa 
} 

/* 

any  uae  of  egatyx  require*  bracketing  by  the  tix  port 
inatruction*  bracketing  the  for  loop  below 
egaqix  taken  from  micro /*y»tema  July 88  20-34 

*/ 

/*  eetting  up  ega  for  write  mode  2  */ 
outportb(0x3oe,5); 
outportb(0x3cf  .2 ) ; 

/*  *et  up  function  reguter  */ 
aam  mov  dx,0x3ce 
aam  mov  al.0x3 
aam  out  dx.al 

/*  function  select  reguter  */ 
aam  me  dx 

/"  react  bit*  3  and  4  */ 
aam  xor  al.al 

/*  set  u  color  atring  pointer,  act  di  s  current  color  */ 
aam  mov  ti .color 
aam  mov  di.  [ti] 

/"  check  for  xor  function  */ 
aam  teal  di.OxlO 
aam  j*  noxor 
/•  act  bita  3  and  4  •/ 
aam  mov  *1,0x18 


noxor 

a am  out  dxul 
/•  bx  a  80  •  y2  */ 
aam  mov  *xjr2 
aam  aal  ax.l 
•am  1*1  ax.l 
hd  aal  ax.l 


i  mov  ba-ax 
i  aal  ax.l 
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SIG  —  NDRAWVLS.C 


asm  >v  ix. 1 

um  add  bi  n 

/*  bx  =  addres  offset,  color  already  in  di  */ 
um  mov  ax.OxaOOO 
aam  mov  ea.ax 

/*  ea  =  ega  ram  starting  addreas  */ 

/*  select  bit  mask  register  */ 
aam  mov  dx ,0x3ce 
aam  mov  al.Oxff 
aam  out  dx.al 
/*  output  bit  maak  */ 
aam  inc  dx 
aam  mov  al.OxSO 
aam  ror  al.cl 
aam  out  dx.al 

/* 

DO  NOT  ASSIGN  VALUES  TO  CX,  SI,  DI  IN  THE  LOOP 

V 

if(cchange) 

for  (;  v2-  >=  yl;  ){  /*  egaqrx(  x.y.colorl:  */ 

aam  mov  ah.e*.[bx] 
aam  mov  ax , di 
aam  mov  ea:[bx],al 
aam  sub  bx,80 
aam  inc  ai 
aam  mov  di.  [si] 

}  /*  for  y2  */ 

else 

for  (;  y2-  >=  yl;  ){  /*  egaqtxf  x.y,  color );  */ 
aam  mov  ah.ea:[hx] 
aam  mov  ax.di 
aam  mov  es:[bx].al 
aam  sub  bx.SO 

}  r  for  yi  V 

o  u  t  portb  ( Ox  3ce ,  5 ) ; 
outportb(0x3cf.0|; 
outportb(  Ox3oe  ,8 ) ; 
outportb(  0x3cf  .Oxffl : 

/*  returning  ega  to  write  mode  0  */ 

}  /*  draw v in  */ 


SIG  —  NSGBLS.H 
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r 

SPEKDtSP:  number  of  vertical  pixels  used  to  dispLsp  power  spectra 

MAXFSIZE:  maximum  fft  site 

HEADWORDS:  site  of  lay  header  in  2-byte  word* 

MAXCOLOR.  the  highest  color  number  used  for  power  ipectral  display 

’BSTCMODE:  video  mode  number  for  the  best  graphics  mode 

7VSCALE:  number  of  vertical  pixels  in  gmode 

7BOTLINE;  last  text  line  number  in  gmode 

STKSZ  the  sise  of  the  FILO  stacks  for  loom /res tore 

HSCALE:  the  number  of  horxcntal  pixels  (same  for  ega,  vga) 

*/ 

#define  SPEKDISP  256 
♦define  MAXFSIZE  256 
# define  HEADWORDS  256 
# define  MAXINT  0X7FFF 
♦define  MAXCOLOR  15 
♦define  VBSTGMODE  0X12 
♦define  VVSCALE  460 
♦define  EBSTGMODE  0X10 
♦define  EVSCALE  350 

♦  define  VBOTLINE  29 
♦define  EBOTLINE  24 

♦  define  STKSZ  6 

♦  define  FDEFAULT  0x2000  /•  Alt-d  */ 

♦define  XNOISE  0x3100  /*  Alt-n  */ 

♦define  FREQDISP  0x2100  /•  Alt-f  */ 

♦  define  DRPLL'S  0x5200  /*  insert  */ 

♦define  L  "’-MINOS  0x5300  /*  delete  */ 

♦aefin  :  Kr  PLUS  0x7700  /*  ctH-home  “/ 

♦define  '  '.MINUS  0x7500  /•  Ctrl- end  •/ 

♦  define  POPGSTACK  0x4900  /•  PgUp  */ 

♦  define  EXITPROG  Oxl lb  /*  Eac  */ 

♦define  PSHGSTACK  0x5100  /*  PgDn  */ 

♦define  WRTFILE  0x7200  /•  Ctrl-PrintSern  */ 

♦define  FFTNCRSE  0x6c00  /•  Alt-F3  */ 

♦define  FFTDCRSE  OxSdOO  /•  Alt-F4  */ 

♦  define  ENTER  13 
♦define  BACKSPACE  8 
♦define  NAM  ESIZE  40 
♦define  CURSCOL  6 
♦define  TTTLCOL  9 
♦define  FCURSCOL  8 
♦define  HSCALE  640 
♦define  MODE-DEFLT  6 

♦  include  <stdio.h> 

♦  include  <io.h> 

♦  include  <fcntl.h> 

♦  include  <stdlib.h> 

♦include  <math_h> 

♦include  <dosh> 

♦  include  <string.h> 


82 


SIG  —  NSIG.PRO 


/'  prototype*  generated  from  function*  m  kay«tufT\fltfunai  c  '/ 

int  fft-»ixe( int  xf); 

int  polaranode(inl  pan); 

float  gain_*et(double  gn); 

void  fldemean(float  'fdata,  int  iaixe); 

void  imean2fp(int  idatQ,  float  rdatQ,  int  itiae); 

void  makeP*pec(float  rdalQ ,  int  Uile); 

void  normPtpec) floet  rdatQ,  float  'avgp,  int  i*ize); 

void  colrPip-ec(float  rdatQ,  char  kolr»Q,  int  iaixe,  float  colacalr ) ; 

int  »cLifTanp(conit  void  *recl,  con*t  void  *rec2); 

float  getacale(int  ifp,  long  fdi*p,  long  ldi*p,  int  xftite,  int  idntQ,  float  rdatQ,  int  *imin.  int  *imar. 
int  xnoiae,  int  h»c*le,  float  avgpQ,  float  *g  jnrtn,  float  ‘aotoiae); 

/'  prototype*  generated  from  function*  in  kay«tufr\dra*hJn.c  '/ 
void  drawHLnfint  y,  int  xl.  int  x2,  char  color); 

/'  prototype*  genoated  from  function*  in  kay»tuff\ndrawvln.c  '/ 
void  draw VLn( int  x,  int  yl.  int  y2,  char  'color,  int  cchange); 

/'  prototype*  genoated  from  function*  in  kny»tufT\aom*ig.c  '/ 
void  ticm*rk*(long  he,  int  vmin,  int  fplot); 

int  linepiot(int  ifl,  long  hntd,  long  hinc,  int  mnval.  int  mxval.  int  fqplotint  xnoi»e); 
void  vcur*e(long  *fln.  long  *Un,  long  fdiip,  long  ldiip.  long  imov,  unaigned  leur*); 
void  hcur*e(int  *ub.  int  *db.  int  imov,  int  upbar); 
void  timcatflong  jcune); 

void  errant  (long  icur*e.  int  row,  int  Utcurie.  int  ip,  long  hnc); 

int  getcur*e*(long  "fvline.  long  “lvline,  int  *ubar,  int  *dbar,  long  *fcur»,  long  'leur*.  tong  hinc, 

int  freqpioi.  int  xnoiae.  int  ifp); 

int  wrtfUe(int  iftle,  long  fdiip.  long  ldiip); 

int  prmchange(un*igned  ivjnt  'tgjnt  *kn); 

float  *cale*et( float  maxpowint  kn); 

void  p*leat(wid); 

/'  prototype*  gensated  from  function*  in  k*y*tu/T\whtunio.c  '/ 
void  errexit(char  *m»g); 

int  idat  read  (int  ip.  int  x*ixe.  long  lakip,  void  *dat); 
long  fopeninp(char  'fname,  int  'ifile); 

void  frxunmax(int  idatQ,  int  nelem.  int  “vniin.  int  *vmax); 

/*  prototype*  generated  from  function*  in  kay*luff\plotutl*.c  */ 
int  readChar(int  page,  char  'atlr); 
int  getkey(); 

v<*d  gotoxy(inl  column,  int  row); 
int  getPage( void): 

void  getxy(int  'column,  int  'row,  int  page); 
void  wm Char( char  ch.  int  color,  int  page); 
int  getMode(int  *ncol»); 
void  cur*orinc(int  « pace*. int  page); 
void  cxlf(in*  page); 

void  wntStringfchar  **tr,  int  color,  int  page); 

void  *etMode(int  mode); 

void  cl* (char  color*); 

vead  readS tnng( char  *«ptr,  int  coimt); 

void  aet palette) in*  palette,  int  color); 

void  aet graphic* (int  gmode); 

void  clriinefint  line); 


SIG  —  PLOTUTLS.C 


83 


#  include  "  nsig.pro" 

#include  <do*.h> 

♦define  BACKSPACE  8 
♦define  ENTER  13 

/* 

READcHAR 

Reeds  4  Cher  ec  ter  from  the  Screen 

7 

int  reed  Cher  (int  pege,  cher  *eur) 

{ 

union  RECS  reg; 
reg.heh  =  8; 
reg.hbh  =  pege; 
inl86  (OnlO.treg.treg); 

*ettr  =  reg.heh; 
return  (reg.h.el); 

} 

/* 

GETKEYO 

uses  ROM-BIOS  eervioe  to  gel  e  keycode  from  buffo-,  or  well  for  keyboard 

V 

int  getkey() 

{ 

union  REGS  reg; 
reg.heh  =  0; 

retumf  int  86(0*  16,t  reg  reg));  /"  turboC  int86()  return*  AX  */ 

} 

r 

GOTOXY() 

Move*  Cursor  to  Specified  x.y  Position  end  Pege  3 

V 

void  golaxy(int  column,  int  row) 

{ 

union  REGS  reg; 
reg.heh  =  3; 
reg.hbh  =  0. 
reg.hdh  =  row; 
reg  hdl  =  column; 
int 86  (Oxl0.treg.treg); 

} 

/* 

GETpAGE() 

Returns  Active  Page  Number 

7 

int  getPege(void) 

{ 

union  REGS  reg; 
reg.heh  =0*0  F . 
int86(  0*10,  t  reg,  treg); 
return  (reg.h.bh); 

} 

r 

GETXY() 

Get*  ths  Current  Cursor  Position 

7 

void  get i y (in t  'column,  int  *row ,  int  pege) 

{ 

union  REGS  reg; 
reg.heh  =  3; 
reg .h bh  »  pege; 
in**6  (OxlO.treg.treg); 

'row  *  reg.h  dh; 

‘colitm  s  reg. hdl; 


} 
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SIG  —  PLOTVTLS  C 


r 

WRTTcHAR 

Write*  *  Character  to  the  Scrcm  with  Specified  Attribute 

*/ 

void  wntChar(dikr  ch.  int  color,  int  page) 

{ 

union  REGS  reg; 
regh-kh  =  9; 
reg. h.  si  =  ch; 
reg.h.bl  =  color; 
reg.h.bh  =  page; 
reg.x.cx  =  1; 
inl86  (OxlO.&reg.ireg); 

} 

/* 

GETmODEO 

Return  Current  Video  Mode 

*/ 

int  getMode(int  *ncob) 

{ 

union  REGS  reg; 
regh-kh  s  OxOF; 
inl£6  (OxlO.fcreg.&reg); 

*ncoU  =  regh-kh; 
return  reg  h.al; 

} 

/* 

CURSORINC 

moves  the  cursor  right  one  space,  wrsps  lines  knd  pkge 

7 

void  cursorinc(spkoes.pkge) 
int  spaces. page; 

{ 

int  col,  row,  width,  length; 

if  ((0x11  ==  (length  =  getMode  (fcwidth)))  ||  (0x12  =  =  length)) 
length  =  30; 

else 

length  =  25; 
get  xy  ( li  col  M  row ,  pkge ) ; 
if  ((col  +=  spkces)  >=  width){ 
col  =  0; 
row+  +  ; 

} 

else  if  (col  <0) 
col  ss  0; 

if  (row  >=  length) 
res*  a  0; 
goto* y( col,  row); 

} 

r 

CRLF 

moves  the  cursor  to  the  left  margin.  down  one  line 

V 

void  crlf(pege) 
int  page; 

{ 

int  col.  row; 

get  xy  (it  col  ,1c  row  .page); 
got  cay  (0  josr  + 1 ); 

} 

/* 

WRTT.TRING 

Writes  k  String  to  the  Screen  with  Specified  Attribute 
can  write  faster  if  cursonnc  functions  are  placed  in  the 


S1G  —  PLOTUTLS.C 
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function  iaitulittiioD,  and  col  and  row  are  updated  within 
the  routine. 

V 

void  wnt StnrLg( char  "*ir,  int  color,  int  page) 

{ 

char  cval;  int  i  =  0; 
while  (cval=»tr[i]){ 
if  C\n'==cval){ 
crlf(page); 

} 

elae{ 

writChar(cval,  color,  page); 
cur>onnc(  1  .page); 

} 

i+  +  ; 

} 

} 

/• 

SETmODEf) 

Set*  Video  Mode 

*/ 

void  *elMode(Lnt  mode) 

{ 

union  REGS  reg; 
reg.h-ah  =  0; 
reghal  =  mode: 
int 86  (OxlO.fcreg.ltreg); 

} 

/* 

CLS()  =  Clear*  the  Screen 

char  colon  =  (B7)bUnk.  (B8-B4)rgb  foreground.  ( B3 ) foreground  intensity, 
(B2-B0)rgb  badtround 
badtround  aeta  the  fill  color  uaed 
colon  haa  a  different  meaning  under  CGA 
taken  from  Graphic*  Programming  mC;  FtT  Steven*  1988 

V 

void  da(di*r  colon) 

{ 

union  REGS  reg; 

in*  column*  .mode. 

mode  =  getMode(XtcoliHnna); 

if  (column*  =  =  80) 

{ 

if  ((mode  ==  0x11)  ||  (mode  ==  0x12)) 
reg  x  dx  —  OxiDAF; 

elae 

reg.x.dx  *  Ox  IMF ; 
regh.bh  =  colon; 

} 

eke 

{ 

pef  j  di  s  0x1838; 

•witch  (color*) 

{ 

earn  1:  reg-hhh  •  0x55; 

break 

caae  2:  reg.hbb  *  OxAA; 

break 

caae  3:  reghhh  *  OxFF; 

break 

default: 

reghhh  =  0; 

break 


} 
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SIG  —  PLOTVTLS.C 


} 

reg.x.ax  =  0x0600; 
reg.xcx  =  0; 
inl86(0xl0,ireg,&reg); 
got<*y(0,0); 

} 

/* 

READeTRING 

Reada  a  ttruvg  from  the  keyboard,  terminated  with  a  \n 

*/ 

#  define  ESCAPE  Oxlb 

void  readS  onng(  char  “»ptr,  int  coul) 

{ 

ini  vidpage; 
vidpage=getPage(); 

whiie(ENTER  !s  (*«ptr++=getkey{)  It  Oxff)) 
i f  (BACKSPACE  ==  *(«ptr-l)){ 
if  (coimt){ 

cunonncf- 1  .vidpage): 

wti  t  Char(  0x20 ,0xf.  v  i  dpage ) ; 

count-; 

•ptr-  =  2; 

} 

eiae 

•ptr-; 

} 

elee  if  (ESCAPE  ==  *(»ptr-l)){ 

*(»plr-l-coum)sr0; 

return; 

} 

elae{ 

wri  t  Cbar(  “( apt  r- 1 )  .Oxf,  vidpage ) ; 
cur*or«nc(  1  .vidpage); 
count -t  +  ; 

} 

cr  If  (vidpage); 

*(»ptr-l )  =  0; 

} 

void  aetpaleUe(int  palette,  int  color) 

{ 

union  REGS  reg; 
reg  h.ah  =  0x10; 
regh-al  =  0; 
reg  h.bh  =  color; 
reg  h.bl  =  palette; 
uuS6  ( Ox  1 0 .  ii  reg .  It  reg ) ; 

} 

r 

SETGRAPHICS 

7 

void  act  graphical  int  gmode){ 
eetMode(gmode); 

»etpai*tt*(1.8)aetpalette(2.1)aetpal*tte(3.33)aetpalette(4.40); 
aetpaktM(5.S)»etpalette(6.4S)setpal«tlc(7.61)aetpalette(8.47); 
aet  pakt  t«(9.S3)  ;aetpalct  te{  10.37)  «ctpalette(  11. 44)»etpaleue(  1236); 
•*tpWette(13.53)aetpaleue<M.J4)aetpal«tte(15,63), 
da((dtar)0x70); 

} 

r 

CLRUNE 

dean  a  line  on  the  diaplagr,  uaing  my  C-bioa  routaoea 

V 

void  clriinefinl  line){ 
it*  i.  P<: 


SIG  —  PLOTUTLS.C 
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PC  =  *etP««e(); 
jottKy(Oime); 

for(i=0u+  +  <80;writCh*r(32,l.pg).cur»orinc<  1  ,pg)); 
gotoByjO  Jjne), 

} 
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SIG  —  WHTIMIOC 


#  in  elude  "aaig-pro” 

#  include  "  negfaU.lt’ 
extern  ovidmode,  vidpage; 

r 

ERREXIT 

7 

void  errexit(char  *mag){ 
fcloaeall) ); 

writStr®g(”error:  "  ,5.  vidpage); 
wri  t  S  t  ring  ( meg ,  5 ,  vid  page ) ; 
delay(JOOO); 

•etMode(  ovidmode); 

-exit(O); 

} 

/* 

IDATREAD 

V 

ini  :datread(u*  ip.  int  xaise,  long  Ukip,  void  *dai){ 
ini  bytearead; 

bytearead  =  read(ip,  dat.  xiixe<<l); 
laeekfip,  lakip<<l .  SEEK-CUR); 
return)  bytearead  >  >  1 ) ; 

} 

/* 

FOPENINP 

V 

long  fopeninp(char  ‘fnarne.  int  *ifile){ 

•tfile  =  open(fnarae,0-RDONLY|0_BINARY); 
while  (-1  ==  •ifile){ 

wrilString("  Source  File  opening  did  not  succeed:  " ,S,  vidpage); 
wn  lSinng(  fnarne,  4,  vidpage  jcrlfl  vidpage); 
writString(”Emer  Ale  name  for  input  file:  ",15. vidpage): 
•traet(fname.O); 
readS  tnngf  fname.O); 
if  (0  ==  -fname){ 

set  Mode(  ovidmode): 

errexjl(" aborted  from  Ale  entry”); 

} 

UNLOAD 

"iftle  =  open(fn«ne,0-RDONLY|O.BINARY); 

} 

Iaeek(*iftle,O.SEEK-END); 

retum(tell(*ifUe));  /•  return  number  of  byte*  •/ 

} 

void  fminma«(mt  idaiQ,  int  nelem,  int  *vmin.  int  *vmax){ 
ini  ival,  i.  imin.  imax,  *iptr; 

for|i=0,  iptr=  idal .  tminximaxx’idat;  ++i<nelem,  ) 
if  ((ivalx’-f+ipir)  >  imax) 
imax  x  ival; 
da*  if  (ival  <  imin) 
imin  x  ival: 

*emax  x  imax; 

aTHUtt  S  imin: 

)/*forjV 
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